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OBJECTIVE  OF  VISUAL  EVOKED  RESPONSE  STUDY 


The  objective  of  the  present  work  was  to  provide  means  to  improve  and 
extend  results  regarding  the  interpretation  of  the  visual  evoked  response 
(VER),  as  reported  in  relevant  literature  and  as  previously  achieved  by  the 
Laser  Effects  Branch,  Radiation  Sciences  Division,  of  the  U.S.  Air  Force 
School  of  Aerospace  Medicine  (USAFSAM).  Our  effort  has  been  undertaken  in 
support  of  the  Air  Force  goal  of  assessing  the  effect  of  visual  stimuli  upon 
the  flier. 


To  assess  the  effects  of  visual  stimuli,  USAFSAM  had  previously  re¬ 
corded  data  from  Macaca  mulatta.  The  data  obtained  were  multi-electrode 
tracings,  covering  primarily  the  area  17  (primary  visual  cortex).  The  ex¬ 
perimental  setup  was  geared  to  obtaining  modulation  transfer  functions  for 
variable  contrast  levels  during  macular  stimulation  (narrow  f ield-of-view) . 

During  the  analysis  of  these  data,  the  USAF  realized  the  inadequacy  of 
conventional  Fourier  Analysis  of  VERs  based  on  the  raw  data;  for  very  long 
recordings  had  to  be  used  to  achieve  a  reliable  estimation  of  the  modulation 
transfer  function.  Similarly,  long  tracings  were  necessary  for  the  infer¬ 
ence  of  the  psychophysical  contrast  level — the  psychophysical  contrast 
threshold,  in  particular. 

Our  objective  was  to  extract  the  relevant  information  from  the  data  by 
alternative  means  and  in  a  more  efficient  manner.  In  other  words,  the  pres¬ 
ent  objective  required  the  tools  of  modern  signal  analysis  which  were  geared 
to  reducing  the  effect  of  various  noise  mechanisms  and  to  enhancing  the 
underlying  signals.  Emphasis  was  to  be  placed  on  schemes  which  would  be 
simple  and  robust.  The  schemes  were  to  be  automated  in  the  form  of 
algorithms  which  would  permit  the  USAF  to  obtain  the  desired  information 
in-house.  [In  this  context,  the  reader  is  referred  to  a  discussion  of  im¬ 
portant  aspects  of  signal  processing  of  VERs  in  the  preceding  work  by 
Gustafson,  Eterno,  and  Jarisch(5).] 

The  material  in  much  of  this  report  is  presented  in  the  following  or- 


(a)  the  background  of  problem; 

(b)  the  highlights  of  results  (e.g.,  improvements); 

(c)  the  processing  of  information;  and 

(d)  the  analysis  which  motivated  the  processing. 

This  format  will  be  of  assistance  to  those  individuals  who  are  only  inter¬ 
ested  in  using  these  algorithms.  These  readers  will  quickly  find  what  they 
can  expect  in  terms  of  improved  signal  extraction  without  working  their  way 
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through  all  the  intricacies  that  led  to  the  particular  processing  method. 
For  the  motivation  behind  these  algorithms,  refer  to  the  technical  details 
following  the  main  results. 

The  one  exception  to  this  format  is  that  we  discuss  the  data  (obtained 
in  previous  research)  from  Macaca  mulatta  before  the  processing  of  human 
subject  data.  Our  reasons  were  that:  we  received  the  data  from  Macaca 
mulatta  first;  and,  these  earlier  data  permitted  a  more  general  discussion 
of  aspects  of  signal  analysis. 


2.  NATURE  OF  'ME  PROBLEM 


The  interest  of  the  Air  Force  was  in  the  study  of:  initially,  animal 
VERs;  and,  later,  human  subject  VERs.  Fo "  obvious  reasons,  the  data  re¬ 
corded  from  human  subjects  were  restricted  to  scalp  electrode  recordings.  By 
contrast,  in  the  animal  preparations,  ep Mural  bipolar  electrode  recordings 
'.with  a  tip  separation,  at  any  given  site,  in  tne  millimeter  range),  were 
possible.  However,  to  emulate  the  human  subject  data,  the  USAF  decided  to 
provide  us,  at  first,  with  epidural  monopolar  recordings  of  Macaca  mulatta. 
These  data  were  expected  to  exhibit  properties  that  might  also  be  encoun¬ 
tered  in  the  human  subject  —  thus  providing,  within  the  given  time  con¬ 
straint,  an  easily  accessible  basis  for  an  evaluation  of  our  processing 
techniques. 

Preliminary  evaluation  of  these  data  by  UHAFSAM  indicated  that,  in  contrast 
to  bipolar  recordings,  they  exhibit  a  complex  relationship  between  the  vis¬ 
ual  stimuli  and  the  estimated  VERs.  The  most  important  difficulties  encoun¬ 
tered  were: 

(a)  high  variability  of  the  VER  under  fixed  experimental  conditions;  and 

(b)  low  signal-to-noise  ratio. 

To  extract  the  desired  information,  we  ha  i  to: 

(a)  separate  the  signal  from  the  noise,  as  well  as  possible;  and 

(b)  reduce  the  variability  through  proper  combination  of  information 

channels . 

The  problems  at  hand  were  further  complicated  by  the  fact  that  well-defined 
quantitative  or  qualitative  models,  for  the  desired  underlying  information 
carrying  signals,  did  not  exist.  The  problem  of  information  extraction  was 
therefore  more  difficult.  Fortunately,  however,  modern  signal  analysis 
addresses  precisely  these  problems  in  a  systematic  fashion. 

For  the  convenience  of  the  reader,  an  example  of  some  data  is  given  and 
is  briefly  discussed.  Data  previously  recorded  from  Macaca  mulatta  and  sup¬ 
plied  by  the  USAF  are  shown  in  Fig.  2.1.  Channels  #1  -  #8  represent  the  volt¬ 
ages  recorded  with  the  arrangement  of  electrodes  as  shown  in  Fig.  2.2; 
channel  #9  merely  represents  the  control  :’or  the  phase  of  the  pattern  rever¬ 
sal.  The  following  features  of  the  tracings  are  striking: 

The  VERs  are  not  discernible  by  the  naked  eye.  The  tracings  contain 
strong  components  of  the  electrocardiogram  (ECG)  in  all  channels. 
Observe  that  these  tracings  coincide,  over  periods  of  several  sec¬ 
onds,  with  a  given  phase  of  the  control  signal  in  channel  #9. 


EEC  SEGMENT  EEC 

Channel  BOUNDARY  —  ■  (?)  SEGMENT 


Example  of  EEG  and  VER  data,  supplied  by  the  U.S.  Air  Force.  Channels  #1-8  correspond 
to  leads  (as  shown  in  Fig.  2.2);  Channel  #2  carries  invisible  VERs;  and  Channel  #9  is 
the  control  signal  of  pattern  reversal. 


All  traces  show  simultaneous  changes  of  similar  patterns  of  the 
electroencephalographic  (£EG)  background  activity. 

The  distinct  EEG  background  patterns  may  belong  to  a  finite  population; 
a  systematic  relationship  also  appears  to  exist  between  the  admin¬ 
istered  drug  level  (Nembutal)  and  the  duration  and  prevalence  of 
certain  of  these  patterns  (not  visible  in  Fig.  2.1). 

A  further  feature  concerning  the  VER,  although  not  visible  in  these 
tracings,  is  also  of  importance.  From  various  studies  we  know  that  the  res¬ 
ponse  to  a  single  stimulus  consists  of  several  wavelets  (4).  Each  of  these 
wavelets  exhibits  different  degrees  of  variability  in  terms  of  amplitude  and 
phase,  typically  increasing  with  added  time  after  the  stimulus.  All  of 
these  features  of  the  data  determined  the  approach  we  had  to  take  and  the 
attention  we  had  to  devote  to  detail. 


3.  APPROACH  TO  REDUCING  THE  VARIABILITY  OF  ESTIMATES  OF  THE  VER 


The  study  of  the  VER  was  particularly  interesting  because  both  Macaca 
mulatta  and  human  data  were  included.  Some  important  differences  exist  in 
the  experimental  conditions  and  in  the  characteristics  of  the  data.  Here, 
we  give  a  brief  overview  of  these  differences.  The  detailed  discussion  of 
the  individual  approaches  and  studies  of  these  data  are  reserved  for  subse¬ 
quent  sections  of  this  report. 


3.1  STIMULUS  CONDITIONS 

The  stimulus  was  provided  by  a  computer-controlled  TV  display.  For  the 
present  experiments  a  colorless  checkerboard  pattern,  of  *1  degrees  squared 
with  4  checks  per  degree,  was  used.  The  reversal  rate  for  stimulating  the 
Macaca  mulatta  had  been  4  Hz;  but, for  the  human  subject  data,  the  reversal 
rate  was  increased  to  15  Hz.  Also,  the  Macaca  mulatta  had  been  stimulated 
monocularly;  but,  for  the  human  subject,  both  eyes  were  used. 

3.2  SUBJECT  PREPARATION 

In  the  earlier  research,  the  Macaca  mulatta  was  sedated  by  Nembutal,  and 
the  pupil  was  dilated  with  atropine.  To  avoid  any  motion  artifacts,  the  eye 
and  the  eyelid  of  the  monkey  were  arrested  with  sutures;  the  eye  was  kept 
moist  with  a  spraying  device.  Conversely,  the  human  volunteers  were  unre¬ 
strained  and  used  binocular  vision. 


3.3  DATA  COLLECTION 

The  data  from  the  Macaca  mulatta  were  recorded  from  5  epidural  elec¬ 
trodes  spaced  horizontally,  10  mm  apart,  to  cover  the  foveal  representa¬ 
tions  on  the  two  hemispheres.  (Note  that  considerable  anatomical  differ¬ 
ences  exist,  between  Macaca  mulatta  and  man,  in  the  mapping  of  the  foveal 
area  into  the  primary  visual  cortex.)  In  the  human  subjects,  3  scalp 
electrodes  were  placed  above  the  (believed)  location  of  the  calcarine  fis¬ 
sure  and  were  spaced  horizontally  30  mm  apart. 

The  data  were  bandpass-filtered  with  a  lower  corner  frequency  of  0.3  Hz 
(6  dB/octave) ;  the  upper  corner  frequency  was  300  Hz  (6  dB/octave)  for  the 
data  from  the  Macaca  mulatta,  and  100  Hz  (12  dB/octave)  for  the  human  sub¬ 
ject  data.  For  the  Macaca  mulatta  the  sampling  frequency  was  512  sam¬ 
ples/sec;  and,  for  the  human  subjects,  256  samples/sec.  (Unfortunately,  the 
upper  corner  frequency  of  the  data  from  the  Macaca  mulatta  was  higher  than 
the  Nyquist  frequency.)  Note  that  the  measurement  of  the  differential  volt¬ 
ages  between  the  electrodes  was  arranged  differently  for  the  two  types  of 
data,  as  explained  in  more  detail  in  section  6.1. 


In  section  4,  we  introduce  the-  data  formerly  obtained  from  the  Macaca 
mulatta.  Then,  in  section  5,  we  present  a  number  of  interesting  findings; 
and,  in  section  6,  we  turn  to  the  modification  of  the  recording  and  analysis 
approach  for  the  human  subjects.  in  section  7,  some  basic  studies  for  these 
subjects  are  given.  Finally,  in  section  8,  we  address  some  ideas  on  improv¬ 
ing  experimental  conditions  for  refining  the  current  modeling  approach. 


If. 


4.  THE  DATA  FROM  MACACA  MULATTA 


Scientific  Systems,  Inc.,  was  provided  with  several  minutes  worth  of 
earlier  data,  obtained  with  and  without  stimulus,  and  sampled  before  and 
after  the  subject  (Macaca  mulatta)  had  been  blinded,  transiently,  by  a 
flash.  Even  though  our  new  data  processing  scheme  has  now  been  tried  with¬ 
out  any  fine  tuning,  considerable  improvement  in  the  estimation  of  visual 
evoked  responses  has  been  achieved.  Furthermore,  we  found  some  interesting 
results  which  suggested  the  usefulness  of  the  combination  of  VER  and  in¬ 
stantaneous  variability  in  assessing  visual  performance. 


4.1  THE  DATABASE 

For  the  development  of  our  basic  software  we  were  provided  with  data 
previously  recorded  by  the  USAF.  These  data  represent  potentials  from 
5  epidural  electrodes,  aligned  horizontally  above  the  visual  cortex.  The 
electrodes  were  spaced  10  mm  apart,  and  potentials  were  recorded  relative 
to  the  left  ear  (Fig.  4.1).  The  data  had  been  preamplified  with  Grass  pre¬ 
amplifiers  (.3  Hz  -  300  Hz),  recorded  on  FM-tape  and,  at  replay,  digitized 
with  512  samples/sec. 

For  the  development  of  our  softwace,  we  concentrated  on  a  dataset  with, 
approximately,  the  following  stimulus  conditions: 


Time 

Conditions 

0 

-  60 

sec 

: 

B3  ank  screen 

60 

-  180 

sec 

• 

• 

Stimulus 

180 

-  420 

sec 

• 

Flash,  followed  by  240  sec  of  stimulus 

The  stimulus  consisted  of  a  4-degree  visual  field  with  a  reversing  checker¬ 
board  pattern.  The  reversal  rate  (rectangular  modulation)  of  the  stimulus 
was  4/sec  at  30%  contrast.  At  T  =  180  sec,  a  .85-Joule  flash  was  dis¬ 
charged,  causing  a  transient  blinding  of  the  subject.  Throughout  the  exper¬ 
iment,  only  the  right  eye  of  the  subject  was  kept  open,  fixed,  and  artifici¬ 
ally  moistened.  Nembutal  was  administered  as  a  sedative. 

We  selected  these  data  from  five  similar  sets,  because  they  provided 
not  only  the  largest  database  but  also  numerous  artifacts  with  which  we 
would  eventually  have  to  cope.  They  appeared  to  provide,  also,  a  more  reli¬ 
able  database  when  compared  with  the  first  set,  as  one  might  expect  various 
gains  settings  to  be  better  adjusted  the  second  time. 


4.2  GRAPHIC  DISPLAY  OF  THE  VER 


An  important  criterion  in  assessing  the  performance  of  data  analysis  is 
the  capability  to  discern  stimulus  conditions.  For  this  purpose,  we  emu¬ 
lated  the  previous  USAFSAM  data  processing  and  compared  it  with  ours. 

Previous  processing  by  USAFSAM  was  based  on  the  amplitudes  of  the  fre¬ 
quencies  of  the  VER,  the  fundamental  at  the  reversal  rate  (here,  4  Hz),  and 
its  harmonics.  However,  we  display  the  data  in  a  slightly  different  way 
than  in  earlier  USAFSAM  studies — rather’  than  plotting  the  amplitude  over 
short  consecutive  data  windows,  per  se,  we  display  the  cumulative  sum  of 
these  amplitudes:  i.e., 


T 

R(T)  =  I  a(t) 


(4.1) 


t  =  0 


(not  power,  since  eq.  4,1  is  more  robust). 

The  idea  here  is  that  a  consistent  change  of  the  amplitude  would  build 
up  an  integral  response,  while  random  changes  would  tend  to  cancel.  How¬ 
ever,  the  amplitudes  will  always  be  positive  quantities; bo*  even  without 
stimulus,  the  foregoing  sum  will  show  some  buildup.  Upon  stimulation,  how¬ 
ever,  one  would  expect  a  faster  buildup  (increase  in  slope).  This  idea  is 
displayed  in  Fig.  4.2. 

Next,  we  proceed  to  the  cumulative  output  of  our  matched  filter  ap¬ 
proach.  The  details  of  this  method  are  explained  in  section  4.4.  At  the 
moment  observe,  in  Fig.  4.3,  the  clear  onset  of  a  stimulus  at  T  =  60  sec, 
compared  with  the  display  in  Fig.  4.2.  This  comparison  gives  the  first  hint 
of  the  way  in  which  we  will  proceed. 

Regarding  Fig.  4.2,  one  may  argue  that  the  onset  of  a  stimulus  is  ob¬ 
scured  because,  visually,  a  change  of  slope  in  a  display  is  hard  to  judge. 

To  cope  with  such  an  argument,  we  produc’od  the  plot  in  Fig.  4.4  (modified 
from  Fig.  4.2),  by  removing  the  average  slope,  which  is  present  without  a 
stimulus.  Although  this  step  appears  to  bring  some  improvement  (additional 
information  has  been  used,  however),  the  integrated  VER  appears  less  clear 
in  Fig.  4.4  than  in  Fig.  4.3. 

Here  observe  also  that,  in  Fig.  4.3,  with  the  onset  of  the  stimulus 
(T  =  60  sec),  the  integral  response  form.’;  a  much  smoother  graph.  The  fact 
that  such  a  change  in  character  is  not  apparent  in  Fig.  4.4  can,  in  itself, 
oe  utilized  to  assess  visual  response. 

In  summary,  by  plotting  the  cumulative  response,  we  demonstrated  the 
superior  performance  of  the  new  matched  filter  approach.  In  section  5,  we 
discuss  the  individual  3teps  of  the  new  ipproach  as  well  as  some  preproc¬ 
essing  of  the  data. 
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Figure  ;K3-  Integral  response  from  the  matched  filter  approach.  Observe 

the  clear  onset  of  a  VER  at  T  -  60  sec.  The  integral  response 
appears  to  become  smoother  once  the  stimulus  is  present.  This 
change  of  the  quality  of  the  response  may  be  of  use  for  assess¬ 
ing  visual  processing. 
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Figure  4.4:  Plots  ji  and  Ji.  Modification  of  the  plots  in  Figire  4.2 
Here,  the  average  slope  of  the  integral  response  during 
a  nonstimulus  period  has  been  removed.  Although  this 
step  enhances  the  onset  of  the  VER,  the  variability  of 
the  response  estimate  remains  high.  Also,  no  change  is 
apparent  in  the  character  of  the  response  curve  (e.g., 
reduced  variability)  . 


4.3  MODELING  APPROACH  FOR  DATA  PROCESSING 


The  philosophy  of  our  data  processing  consists,  in  general,  of  three 
parts: 

(a)  removal  of  artifacts  (modeling  of  artifacts); 

(b)  transformation  of  data  (modeling  of  EEG) ;  and 

(c)  extraction  of  signal  (modeling  of  VER). 


•4.3.1  Removal  of  60  Hz  and  Harmonics 

Power-line  noise  (60  Hz)  and  its  harmonics  impose  some  interferences  in 
CEG/VER  recording.  For  better  analyses  of  the  VER  to  be  performed,  these 
artifacts  have  to  be  removed  from  the  raw  data.  A  robust  adaptive  filter 
has  been  developed  to  achieve  this  purpose.  The  mechanics  of  the  filter  are 
quite  straightforward.  First,  a  trigonometric  identity  is  employed: 

sin(0+6)  =  2sin0cos6  -  sin(0-6)  .  (4.2) 


It  represents  the  phase  change  between  digital  samples  for  a  sinusoidal 
signal.  Equation  4.2  can  be  used  to  calculate  the  sinusoidal  amplitude  at 
the  next  sample  time  (phase  equal  9+6)  in  terms  of  the  amplitude  at  the 
previous  two  samples  (phases  of  6  and  0-6).  Multiplication  of  (4.2)  by  the 
amplitude  magnitude  results  in: 

a*+1  -  2aiCos6  -  a^-i  (4.3) 

where  represents  the  time  series  of  alternating  current  (a.c.)  interfer¬ 
ence  amplitudes.  If  we  denote  the  digitized  source  signal  sequence  by  > 
then  the  difference,  x^-a^,  should  be  a  null  apart  from  a  possible  constant 
offset  (in  our  case,  this  offset  is  EEG/VFR  background  at  sample  i).  The 
previous  difference,  xj-i  -  aj_i ,  is  the  best  available  estimate  of  this 
offset.  In  equation  4.4: 

f  =  (xi-aj)  -  (xi-i  -  ai-i )  ,  -  (4.4) 

The  quantity  f  will  be  null  if  the  interference  amplitude  is  properly  esti¬ 
mated  by  a ], .  Thus,  f  should  be  used  as  a  positive  feedback  element.  Based 
on  equation  4.3,  the  logic  followed  is: 

if  f  >A ,  increase  aj  by  a  small  fixed  amount, 

if  f<-&,  decrease  a^  by  a  small  fixed  amount, 

where  A  is  a  small  positive  number. 
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This  logic  is  implemented  in  l.iG/VER  analysis,  with  the  results  shewn 
in  Fig.  4.5.  The  result  is  quite  >  3bust  in  the  sense  that,  if  a^  is  overes¬ 
timated  or  underestimated,  f  will  ‘.end  to  be  negative  or  positive  ,  thus 
causing  to  be  decremented  or  incremented  in  a  cautious  fashion.  The 
feedback  loop  is  stable.  In  addition,  the  adaption  of  the  filter  to  chang¬ 
ing  condition  is  reasonably  rapid. 


4.3.2  Effect  of  Removal  of  ECG 

In  our  data  (FILE  2),  the  ECG  Ls  a  prominent  feature  on  all  five  chan¬ 
nels.  The  spectrum  of  the  nonperiodic  ECG  extends  from  about  2  Hz  well  into 
the  80-Hz  region.  From  some  spectral  analysis,  we  concluded  that  energy  in 
the  20*  to  40-Hz  region  should  provide  good  information  (high  signal-to-noise) 
for  the  generation  of  ECG  fiducial  points.  Thus,  we  decided  to  look  at 
bursts  of  energy  around  30  Hz.  We  use  the  scheme  shown  in  Fig.  4.6  for  the 
detection  of  R-waves.  In  order  to  avoid  confusion  with  P-waves,  we  check 
over  a  200-ms  window  for  the  presence  of  other  bursts  of  energy;  if  none 
occur,  the  detected  local  maximum  i.3  associated  with  a  data  window  from 
which  the  average  ECG  is  computed  (Fig. 4. 7). 

A  second  pass  through  the  data  subtracts  the  average  ECG  whenever  its 
presence  is  detected.  In  Fig.  4.8,  one  can  see  the  effect  of  this  process: 
after  an  initial  startup  period,  the  ECG  and  the  60  Hz,  plus  its  harmonics, 
are  removed  from  the  data. 

The  power  at  the  fundamental  frequency  and  low-order  harmonics  is  com¬ 
parable  or  even  larger  than  the  signal  of  interest — the  VER  (Fig.  4.9).  Ob¬ 
serve  that  a  slight  change  in  heart  rate  could  result  in  a  strong 
interference  of  the  ECG  with  the  VEK.  For  this  reason,  the  ECG  had  to  be 
removed.  A  typical  result  of  such  removal  is  shown  in  Fig.  4.10.  Note  the 
drastic  change  of  the  amplitude  of  the  fundamental  frequency  (2  Hz)  and  the 
first  harmonic  (4  Hz),  next  to  the  fundamental  reversal  rate  of  the  VER. 

Thus,  our  scheme  appears  to  provide  a  good  safeguard  against  any  interfer¬ 
ence  of  the  ECG  with  the  VER. 


4.3.3  Method  of  Construction  of  an  ECG  Template 

The  construction  of  the  template  for  the  ECG  requires  some  pre¬ 
processing  of  data,  and  is  based  on  an  iterative  refinement  of  esti¬ 
mation.  Thus,  we  pass  the  raw  data  first  through  a  high-pass  filter  to 
extract  the  high-frequency  component,  of  the  ECG.  Any  clipping  of  the  data 
due  to  the  analog-to-digital  (A/D)  converter  is  currently  detected  in  the 
high-pa33  filtered  version  of  the  data;  and  the  data  corresponding  to  clip¬ 
ping  are  excluded  from  further  processing.  Following  this  signal  condition¬ 
ing  .preliminary  fiducial  points  for  the  ECG  are  computed  by  use  of  the 
the  accumulated  velocity,  v. 


Figure  4.6.  Scheme  for  the  determination  of  the  ECG  template.  This  scheme 
does  not  depend  on  the  shape  of  the  F.CG,  since  powe r  in  the 
neighborhood  of  30  Hz  is  observed.  The  200-ms  window  serves 
the  rejection  of  loca.'  maxima  induced  by  P-waves;  thus,  trig¬ 
ger  points  are  associated  with  R-waves. 


Figure  4.7:  Plots  _a  and  Jj,.  Example  of  the  average  ECG  in  channel  #  1 
(  ■  template)  and  standard  deviation  of  ECG  (without  re¬ 
moval  of  baseline).  From  the  shape  of  these  two  plots, 
one  can  estimate  a  timing  uncertainty  for  the  detection 
of  the  R-wave  in  the  order  of  3-ms  standard  deviation. 
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igure  U.8.  The  effect  of  the  removal  of  the  60-Hz  noise  and  the  ECG.  The 
graphs  show  the  processed  data  before  and  after  the  activation 
of  the  60  Hz  and  ECG  removal. 
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Plots  ja  and  Jj.  Comp.irr.son  of  spectral  components  of  the 
VER  and  the  ECG .  (Observe  the  dominance  of  ECG  spectral 
components.) 
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Figure  4.101  Plots  and  b^.  Comparison  of  spectral  components  of 
VER  and  ECG  after  removal  of  the  average  ECG,  (Note 
the  suppression  of  the  ECG  fundamental  and  first  har¬ 
monics  components.) 


Thus,  we  use  the  five  velocities  oi  the  EGG  in  the  channels  by 
computing: 


I  w  i  •  v  i 
i  =  l  —  5 


(4.6) 


For  these  data  the  relevant  weight,-;  ,w,  are  chosen  by  hand,  on  the  basis  of 
subjective  confidence  in  the  clari'.y  of  the  ECG-wavelets. 

Next,  a  preliminary  template  3  computed  for  each  channel  by  use  of 
these  preliminary  fiducial  points.  The  templates  are  moved  over  the  data 
and  adjusted  with  respect  to  their  shape  and  fiducial  point,  as  shown  in 
Fig.  ii.1l. 
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Figure  il.11.  Estimation  of  fiducial  point  and  template. 


4.3.4  Clipping  Of  Data 

All  programs  disregard  data  points  associated  with  clipping  and  (mainly 
at  the  QRS-complex)  adjust  the  weight  of  estimates  accordingly.  For  our 
data,  the  estimates  derived  in  this  v/ay  (by  regarding  clipped  data  as  miss¬ 
ing)  are  biased  estimates;  but  the  bias  tends  to  be  small  compared  to  the 
sample  variability  and  does  not,  at  this  time,  justify  any  correction. 

4.3.5  Modeling  the  EE G-  and  Segmenting 

Visual  examination  of  data  which  are  free  of  60  Hz  and  the  ECG  indi¬ 
cates  that  a  second  autoregressive  model  might  allow  a  good  description  of 
segments  of  the  EEG.  From  our  preliminary  studies,  we  found  that  a 


uonstat ionary  first-order  model,  of  the  type: 


n  =  yt-i+  Gt 

E  [et]  =  0 

2 

var[et]  =  a  for  segment  #  k 


(4.6) 


provides  a  statistically  good  approximation  (that  is,  on  95%  confidence 
level,  e.g.,  only  5%  of  all  models  are  not  well  described)  for  a  few  seconds 
worth  of  data.  The  distribution  of  the  autoregressive  coefficients  and 
«2  is  given  in  Fig.  4.12.  A  clustering  around  ot^  =  1,  c<2  =  0  (equivalent 
to  eq.  4.6),  is  apparent. 


We  note  that  these  coefficients  were  estimated  by  regarding  as  missing 
data  the  clipped  data  and  data  near  R-waves  of  the  ECG.  We  note  also  that 
the  data  displayed  correspond  to  a  sampling  rate  of  256  samples/sec.  The 
variance  of  the  EEG  appears  to  form  two  groups:  one  with  very  large  vari¬ 
ance,  and  one  with  little  variance  (Fig.  4.12).  The  estimate  of  the  vari¬ 
ance  within  each  segment  of  the  EEG  is  anticipated  to  permit,  later,  the 
proper  weighing  of  the  estimated  VER  in  order  to  minimize  variability  of 
response  estimates. 


4.3.6  Generating  a  VER  Template 

The  generation  of  a  template  for  the  VER  is  based,  at  this  time,  on  av¬ 
eraging  in  synchrony  with  the  reference  channel.  When  clipping  occurs,  cor¬ 
responding  data  sections  are  again  treated  as  missing  data  (Fig.  4.13  ). 


4.3.7  Smoothing  the  VER  Template 

From  simple  models  of  the  VER,  we  know  that  only  integer  multiple  fre¬ 
quencies  of  the  reversal  rate  (RR)  can  he  contained  in  the  VER.  For  this 
reason,  we  select  from  the  averaged  VER  only  those  characteristic  frequen¬ 
cies  to  generate  a  smoothed  VER  template.  As  a  processing  option,  the 
number  of  harmonics  can  be  selected.  So  far,  we  select  an  RR  of  4  Hz  and 
harmonics  up  to  16  Hz.  Thus,  in  this  case,  the  VER  is  represented  by 
8  parameters — 2  for  each  characteristic  frequency  (real  and  imaginary  parts). 
The  result  of  this  smoothing  is  shown  in  Fig.  4.14.  This  smoothed  template 
is  used  as  the  reference  for  the  matched  filter. 


4.4  THE  MATCHED  FILTER  APPROACH  FOR  MAC AC A  MULATTA 

A  well-known  fact  is  that  the  matched  filter  is  optimal  in  the  least 
squares  sense  for  estimating  the  (unknown)  amplitude  of  a  waveform  with 
known  shape  embedded  in  white  noise.  Shown  in  Fig.  4.15  is  the  data  proc¬ 
essing,  using  this  matched  filter,  substituting  the  estimated  and  smoothed 
VER  for  the  "known"  waveform. 
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.13:  Plots  a_  and  b.  A  template  of  the  VER  obtained  by 

synchronous  averaging  using  preprocessed  data  (free  of  I 

F.CG  and  60  Hz) .  Clipping  of  data  has  been  taken  into 

account. 
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Figure  A. 14.  The  smoothed  VER  template,  after  removal  of  higher  harmonics 
and  of  all  frequencies  other  than  a  multiple  of  the  reversal 
rate. 
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Figure  4.15.  The  matched  filter  approach  applied  to  the  VER  signal  for 
estimation  of  its  (Integrated)  amplitude. 


To  generate  the  necessary  white  background  noise,  we  transform  our  data 
to  the  innovation  sequence  generated  from  estimated  autoregressive  models. 
From  preliminary  studies  of  segmenting  we  simply  choose:  ai  =  .99,  ap  =  0; 
or,  for  comparison,  a;  =  02  =  0.  From  visual  examination  of  results  we 
cannot,  at  the  present  time,  see  a  significant  change  of  the  estimated 
strength  of  the  VER  (aside  from  differences  in  some  fine  detail  associated 
chiefly  with  the  finer  resolution  display  of  .5  sec).  Compare  the  result  of 
the  integral  response  of  Fig.  4.3  with  that  of  Fig.  4.16.  Note,  however, 
that  the  two  templates  which  were  matchec  against  the  ( transformed)  data 
look  quite  different  in  these  two  cases  (Fig.  4.17). 

Two  observations,  regarding  outputs  of  the  integral  response  of  the 
matched  filter,  deserve  attention: 

First,  one  can  observe  a  reduced  variability  of  the  output  when  the 
subject  was  stimulated.  This  observation  suggests  combining  the  variance 
and  the  amplitude  estimate  for  assessing  stimulus  condition.  Some  more 
spectral  estimates,  for  comparing  the  EEC  background  (one  of  the  components 
introducing  variability  into  estimates)  with  and  without  stimulus,  will  be 
useful  in  clarifying  our  observations  further. 

A  second  observation  (Fig.  4.18)  deserves  further  study — the  recovery 
of  the  VER,  some  10  sec  after  the  flash  at  T  =  180  sec,  appears  to  be  rather 
abrupt.  Not  much  of  a  "slow"  recovery  is  detectable.  This  behavior  is 
reminiscent  of  the  recovery  of  an  amplifier  with  feedback  gain  control  after 
overdriving  it  with  an  impulse. 


4.5  CONCLUSIONS  ABOUT  DATA  FROM  MACACA  MULATTA 

Removal  of  the  EEG  background  from  the  signal  channel  by  use  of  neigh¬ 
boring  channels  (as  suggested  in  our  original  proposal)  turned  out  not  to  be 
useful,  because  the  VER  is  contained  in  all  channels.  Hence,  by  use  of  the 
originally  proposed  regression  model,  one  would  not  only  eliminate  EEG  back¬ 
ground  but  also — to  an  unpredictable  degree — the  VER.  Thus,  we  proceeded 
with  a  separate  processing  of  individual  channels.  Furthermore,  we  have 
reason  to  believe  that  fewer  electrodes  collect  most  of  the  relevant  infor¬ 
mation;  e.g.,  using  only  the  outputs  from  a  pair  of  channels  from  the  left 
and  right  hemispheres'  may  be  suf f’icient*  The  combination  will  depend  on  the 
correlation  structure  of  the  outputs  of  our  matched  filters  (one  for  each 
hemisphere) . 
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Figure  4.16.  Result  of  matched  filler  with  ai  =  .99,  ap  =  0.  The  integral 
response  is  plotted  with  .5  sec  resolution. 
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Figure  4,17:  Plots  a^  and  b^.  Comparison  of  transformed  template  with 
original  template. 
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5.  ANALYSIS  OF  RAW  DATA  FROM  MACACA  MULATTA 


Some  Basic  Properties 


Before  we  began  designing  our  data  processsing  scheme,  we  analyzed  our 
data  by  spectral  analysis,  by  auto-  and  cross-correlation,  and  by  spectral 
methods.  Even  though  our  eventual  data  processing  scheme  is  rather  a  time 
domain  estimation  technique,  we  found  tho  spectral  analysis  to  be  more  re¬ 
vealing  for  assessing  gross  properties  of  the  data,  mainly  because: 

a.  The  VER  is  a  periodic  signal. 

b.  The  noise  introduced  by  the  EEG  is  a  near-periodic  process. 

c.  The  60-Hz  noise  and  its  harmonics  provide  relatively  pure  spectral 

lines. 

Our  observations  concerning  these  processes  are  presented  in  the  following 
material . 

5.1  COVARIANCE  STRUCTURE 

At  least  three  components  of  the  signal  are  of  interest,  and  these  are 
shown  in  Figures  5.1  -  5.3: 

a.  a  slow  component  with  a  period  cf  some  3  sec; 

b.  a  strong  serial  dependency  of  tie  ECGs;  and 

c.  a  strong  60-Hz  signal  (and,  presumably,  harmonics). 

The  slow  component  could  be  associated  with  the  .3-Hz  filter  cutoff 
frequency  prior  to  A/D  conversion,  or  could  be  physiological  (e.g.,  related 
to  respiration) . 

Next,  the  ECG  serial  dependency  is  interesting,  as  it  suggests  the 
possibility  of  a  better  way  to  remove  the  ECG  by  modeling  its  serial 
dependency  with  a  simple  model,  perhaps  by  a  low-order  autoregressive 
model . a 

Finally,  observe  the  strong  60-Hz  signal  (and  its  harmonics).  Such 
strong  signals  have  to  be  removed. 

5.2  PROPERTIES  OF  THE  VER  AND  THE  ECG 

As  a  first  step,  we  tried  to  assess  the  properties  of  the  VER.  For 
this  purpose,  we  examined  the  coherence  between  stimulus,  x,  and  the  recorded 
brain  potential,  y  (Fig.  5.4).  That  is,  we  estimated: 


a  Note:  An  example  of  the  underlying  ECG  (in  terms  of  an  estimated 
template)  is  shown  in  Figure  5.6. 
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Up  to  the  third  harmonic,  (as  can  be  seen)  ,  our  data  show  significant 
coherence  (channel  2  of  File  2).  'hus,  when  smoothing  the  VER,  we  allow  up 
to  this  harmonic  in  our  template. 

Next,  we  examined  the  coherence  of  the  EEG  between  the  electrodes. 

Our  result  can  be  summarized  easil/:  Within  each  hemisphere,  high  coherence 
is  observed  up  to  high  frequencies;  e.g.,  40  Hz,  C(f)>.5.  In  contrast, 
the  coherence  between  the  left  and  right  hemispheres  is  limited  to  low  fre¬ 
quencies;  e.g.,  above  8  Hz  C(f) <.2  (Fig.  5.5). 

These  observations  suggest  that  most  information  about  the  VER  is 
carried  in  any  pair  of  electrodes  of  opposite  hemispheres.  Adding  any  other 
channel  will  not  provide  much  new  information.  Thus,  as  far  as  data 
acquisition  and  analysis  are  concerned,  one  may  concentrate  on  only  a  pair 
of  channels  from  opposite  hemispheres  —  without  losing  much  information 
because  of  disregarding  other  channels. 


5.3  THE  SPECTRAL  PROPERTIES  OF  THE  ECG 

The  electrocardiogram  consists  of  several  wavelets,  principally  the 
P-wave,  the  QRS-complex,  and  the  T-wave.  All  of  these  components  are  easily 
discernible,  especially  from  the  average  of  waveforms  (Fig.  5.6).  Most  of 
the  energy  is  contained  in  the  R-wave.  This  short  impulse-shaped  wave 
generates  harmonics  of  the  heart  rate  well  above  60  Hz. 

Because  of  this  wide  spectrum  (Fig.  5.5),  some  freedom  exists  in  se¬ 
lecting  a  frequency  window  for  the  detection  of  the  ECG,  which  is  free  of 
harmonics  of  the  VER  and  of  60-Hz  noise.  From  a  rough  analysis  of  the  aver¬ 
age  ECG  and  its  standard  deviation  (Fig.  5.6),  we  estimate  the  uncertainty 
in  the  time  of  detection  of  the  ECG  to  be  about  3  ms,  which  is  slightly  more 
than  the  sampling  rate. 


5.4  THE  60-Hz  NOISE 

Our  data  had  been  sampled  at  512  samples/sec.  This  sampling  rate  al¬ 
lowed  us  to  detect  a  60-Hz  component  with  its  harmonics  i’*>  to  300  Hz  (which 
appears  aliased  at  212  Hz).  Thus,  we  wish  to  emphasize  c.^»t,  in  future 
digitization,  the  high-frequency  cutoff  of  the  filters  has  to  be  set  to  much 
lower  frequencies  and  to  a  steepei-  rolloff. 
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Figure  5.6:  Graphs  &  and  b_.  Average  ECG  (AECG)  and  standard  deviation 
of  the  ECG  (SECG)  in  channel  2. 


Another  observation  was  of  interest:  the  60-Hz  noise  has  sidebands  at 
60 1  4  Hz,  coherent  with  the  VER.  These  sidebands  could  be  due  to  the 
nonlinearities  of  the  recording  equipment,  or  to  some  variation  of  the 
impedance  of  the  visual  cortex  in  synchrony  to  the  VER.  At  the  time,  how¬ 
ever,  we  could  not  pursue  this  lead  any  further.  One  way  to  discriminate 
between  the  two  possibilities  is  to  determine  if  the  strong  ECG  also 
modulates  the  60-Hz  noise.  If  so,  the  recording  equipment  must  be  held 
responsible  for  these  modulations. 

If  the  ECG  does  not  modulate  the  60-Hz  noise,  then  the  recording  equip¬ 
ment  behaves  linearly  and  cannot  be  responsible  for  modulation  of  the  60-Hz 
noise;  possibly,  impedance  variations  of  the  visual  cortex  are  detectable. 

5.5  CONCLUSION 

The  potential  for  a  number  of  techniques  to  improve  the  signal-to-noise 
(S/N)  ration  of  the  VER,  versus  other  electrical  activities,  became  apparent 
from  these  studies.  Three  aspects  of  thene  data  made  the  analysis  more 
difficult,  however,  than  was  anticipated  in  the  original  proposal: 

a.  The  VER  was  contained  in  all  channels,  not  only  in  those  placed 

above  the  foveal  projections  in  area  17. 

b.  Clippings,  mainly  associated  with  the  occurrence  of  “he  QRS-complex 

of  the  ECG,  occurred  frequently  in  some  of  the  channels  and  com¬ 
plicated  removal  of  the  ECG. 

c.  A  powerful  60-Hz  component  and  its  harmonics  were  encountered  in  all 

channels,  and  required  appropriate  prefiltering. 

Despite  these  complications,  significant  enhancement  of  the  extraction 
of  the  VER  was  possible.  The  analysis  of  these  USAFSAM  data  provided, 
therefore,  a  good  basis  on  which  to  begin  analyzing  the  human  subject  data. 
Because  the  human  subject  data  were  yet  to  be  recorded,  the  results  of  the 
foregoing  analysis  and  the  increased  understanding  of  the  nature  of  the 
problem  could  be  used  for  improving  the  recording  technique.  In  fact,  the 
modified  approach  resulted  in  data  of  much  better  quality.  In  the  following 
report  sections,  we  address  these  human  data,  and  demonstrate  the  signifi¬ 
cant  improvements  that  are  possible. 


6.  HI  MAM  SUBJECTS 


In  this  phase  of  the  research,  we  used  human-scalp  electrode  tracings. 

For  the  recording  of  these  data  (by  the  present  principal  investigators  at 
the  USAFSAM)  some  important  results  from  the  previous  epidural  recordings  on 
Macaca  mulatta  could  be  used  in  a  beneficial  manner. 

The  approach  used  for  human  suoject  data  analysis  follows,  in  general, 
the  outline  of  our  proposal.  In  summary,  a  rather  drastic  reduction  of 
variability  in  response  estimates  has  been  attained,  especially  in  the  sub¬ 
ject  with  high  variability.  Impressive  results  have  been  achieved  despite 
the  fact  that  the  data  exhibited  a  jehavior  somewhat  different  not  only  from 
the  data  available  to  us  during  the  development  of  the  proposal,  but  also 
from  the  data  recorded  on  the  Macaca  mulatta. 

The  present  report  compares  the  current  results  of  Scientific  Systems, 
Inc.,  with  the  earlier  USAFSAM  results.  The  effectiveness  of  the  individual 
layers  in  our  approach  is  demonstrated.  Then  we  turn  to  the  details  of  the 
current  processing  method. 

Finally,  we  discuss  some  of  the  interesting  properties  of  the  data. 

These  properties  —  which  have  been  unraveled,  in  part,  by  separate  statis¬ 
tical  analysis  of  the  data  —  show  that  the  concept  of  variability  has  to 
be  viewed  in  a  new  light.  For  example,  variability  in  itself  is  a  predictor  of 
stimulus  condition  and  has  a  possible  relation  to  hemispheric  differences 
between  the  sexes  (1,  2,  and  3).  We  now  begin  our  comparison  of  recording 
techniques. 


6.1  RECORDING  TECHNIQUES 

The  first  step  in  recording  is  the  placement  of  electrodes  and  the 
specification  of  the  potential  differentials  which  are  to  be  measured.  In 
the  past,  USAFSAM  used  the  arrangement  shown  in  Figure  6.1.  The  first  set 
of  data  we  had  analyzed,  recorded  epidurally  from  Macaca  mulatta,  was  based 
on  that  arrangement.  Our  analysis  of  that  data  revealed  ( in  contrast  to 
earlier  information  given  to  us)  the  presence  of  a  VER  on  all  electrodes. 
Furthermore,  we  found  the  responses  to  be  nearly  identical  In  the  corre¬ 
sponding  channels,  much  like  a  common  mode.  This  situation  results  in  a 
number  of  limitations  when  analyzing  and  interpreting  VERs. 

To  combat  such  limitations,  we  set  out  to  record  in  an  arrangement  from 
which  we  could  anticipate  good  common  mode  rejection  and  thus  facilitate 
analysis  of  data.  The  new  arrangement  (Fig.  6.2)--which  assigns  spatial 
differentials,  relative  to  the  center  electrode,  to  individual  channels 
(rather  than  relative  to  the  ear) — proved  also  to  be  beneficial  in  regard  to 
avoiding  drift  phenomena  that  lead  to  saturation  of  the  preamplifier  and 
FM-tape.  An  example  of  data  recorded  with  this  arrangement  is  shown  in 
Figure  6.3. 


A  second  refinement  of  our  recording  technique  (also  suggested  by  the 
analysis  of  earlier  USAFSAM  data)  was  the  high-pass  prefiltering  of  the  data 
in  the  preamplifiers.  These  prefiltered  fata  formed  a  second  set  of  data 
(linearly  dependent  on  the  unfiltered  data).  Quite  obvious  during  the  re¬ 
cording  was  the  fact  that  this  technique  avoided  any  risk  of  saturating  the 
recording  media  due  to  electrode  drift.  Later,  during  analysis  of  these 
data,  we  confirmed  an  additional  simplification  with  regard  to  the  process¬ 
ing  of  these  data  (Figs. 6.4  and  6.5). 

The  new  human  subject  data  were  recorded  on  FM-tape,  along  with  the 
1 5-Hz  reversal-rate  reference  signal  and  n  additional  information  channel 
for  the  EMC  associated  with  eye  blinking.  The  stimuli  consisted  of  10  repe¬ 
titions  of  the  following  experiment:  The  subject  was  seated  in  a  dark, 
electromagnetically  shielded  room  in  front  of  a  computer-controlled  TV  dis¬ 
play.  At  first,  60  sec  of  blank  screen  was  displayed,  followed  by  80  sec  of 
a  checkerboard  pattern  at  30-percent  contrast  and  15-Hz  reversal  rate.  The 
visual  field  was  about  4  degrees,  squared,  with  about  4  checks/degree.  With 
these  data  we  can  demonstrate  a  drastic  reduction  of  "variability"  as  com¬ 
pared  with  the  earlier  USAFSAM  analysis. 


6.2  COMPARISON  OF  VARIABILITY  OF  VER  ESTIMATES 

Our  original  idea  of  reducing  variability  in  response  estimates  is 
based  on  extracting  the  information  from  the  relevant  information  channel, 
more  efficiently  than  by  the  original  FFT  method  used  by  USAFSAM.  Further¬ 
more,  "neighboring"  channels  should  be  used  for  reducing  the  effect  of  the 
REG  background  activity.  As  mentioned  before,  however,  all  channels  tend  to 
contain  VERs.  Hence,  information  contained  in  these  channels  is  also  a 
mixture  of  VERs  and  EEG  background  activity. 


To  deal  with  this  situation,  we  have  modified  our  approach  slightly. 

We  treat  all  the  individual  channels  with  our  matched  filter  before  combin¬ 
ing  their  outputs  by  a  regression  technique.  [Note:  Our  original  idea  was, 
first,  to  use  regression  to  remove  background  EEG  activity  from  the  informa¬ 
tion  channel  (carrying  the  VER);  and  secor.d,  to  apply  a  matched  filter.] 


On  the  basis  of  our  current  understanding  cf  the  VER,  it  can  easily  be 
shown  to  be  the  current  approach  suboptirrf.  1 — typical  for  such  hierarchical 
schemes.  (For  a  summary  of  our  data  processing,  refer  to  Figure  6.6.) 


Despite  the  suboptimal  quality  of  the  current  scheme,  we  did  achieve 
drastic  improvements.  To  demonstrate  tne  effect  of  the  individual  process¬ 
ing  layers  in  the  hierarchical  approach,  we  begin  with  a  short  outline  of 
the  FFT  method  by  USAFSAM,  followed  initially  by  the  effect  of  the  matched 
filter  and,  finally,  by  the  effect  of  integration  of  the  outputs  of  the 
matched  filters. 


and  right  hemispheres  can  easily 
be  distinguished.  Note  that, 
here,  scalp  electrodes  were  applied 


reference 


presence  of  15-Hz  reversal-rate  stimulus 


center  electrode 


channel  9  4  5  6 


Figure  6.4.  Specification  of  channels.  HP  indicates  a  high-pass 
filter  which  served  preprocessing  of  the  VER.  Thus, 
compared  with  previous  USAFSAM  recording,  data  are  dif¬ 
ferentiated  in  space  and  time. 
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Figure  6.5.  Example  of  high-pass  filtered  signal  from  center  electrode. 

Obviously,  electrode  drift  is  no  problem,  and  saturation  of 
amplifiers  can  easily  be  avoided  while  retaining  the  infor¬ 
mation  in  the  signal. 
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Figure  6.6.  Summary  of  TiR  processing 


6.2.1  The  FFT  Method 


The  FFT  method  finds  the  respoise  amplitude  of  the  VER  at  the  reversal 
rate — in  this  particular  instance,  it  the  1  5-Hz  reversal  frequency.  The 
actual  computation  is  simplified  by  taking  short  segments  of  data  (in  our 
case,  of  1-sec  duration),  and  averagirg  a  number  of  them  in  synchrony  with 
the  reference  signal.  For  our  comp  i.-’ison  studies,  40  sec  of  data  were  thus 
averaged  by  averaging  40  segments  o'  cata.  Following  the  averaging  process, 
the  15 -Hz  amplitude  was  computed. 

The  result  for  the  female  subject  is  shown  in  Figure  6.7  (Plot  a),  and 
is  labeled:  FFTDJA1 .  The  display  snows  the  response  estimates  for  the  10 
experiments — first,  giving  the  10  estimates  obtained  from  the  no-stimulus 
situation;  and,  second,  the  10  stimulus-present  situations.  The  display 
shows  the  results  for  the  center  electrode  (channel  1),  which  is  believed 
to  give  some  average  response  from  the  two  hemispheres.  To  facilitate  com¬ 
parison  of  the  FFT  method  with  our  method,  we  transformed  the  raw  FFT  out¬ 
puts  by  removing  the  bias  of  the  estimate  (which  is  the  non-zero  estimate 
present  when  no  stimulus  is  applied!;  and  then  we  scaled  the  display  to  give 
an  average  unit  response  when  the  stimulus  is  present.  When  the  necessary 
scaling  coefficients  are  computed  in  this  way,  one  obtains  (for  the  sequence 
of  10  successive  experiments)  the  outputs  shown  in  Appendixes  1.1  -  1  . *1 .  In 
these  Appendixes,  observe  the  increase  of  the  standard  deviation  SIG-Y  of 
the  estimate  Y-HAT  when  a  stimulus  s  present!  The  sequence  of  10 
sequential  estimates  listed  in  Appendix  1  are  the  estimates  based  on  the  FFT 
data  up  to  (inclusive)  the  particular  transformed  FFT  output,  Y-HAT.  y-BAR 
is  the  average  of  the  previous  plus  current,  Y-HAT. 

For  accuracy  we  point  out  that  at  execution  of  the  FFT  method — due  to 
an  oversight — 5  percent  more  data  wore  actually  averaged  than  intended. 

This  percentage  corresponds  to  42  soo  of  raw  data  for  each  displayed  data 
point:  one  extra  second  of  raw  data  had  been  added  at  the  beginning  and  at 
the  end  of  the  data  section. 

The  analysis  of  the  individual  lata  sections  for  each  of  the  10  experi¬ 
ments  began  at  about  10  sec  after  the  onset  of  the  experiment,  when  no 
stimulus  was  present;  or,  alternatively,  at  about  80  sec,  when  a  stimulus  was 
present. 

The  equivalent  results  of  the  analysis  for  the  male  subject  are  shown 
in  Figure  6.8:  Plot  a,  and  are  labe  ed:  FFTKJA1 .  Observe  the  high  variabil¬ 
ity  of  the  response  estimate  for  th  a  individual  when  compared  with  that  for 
the  female  subject. 

Having  clarified  our  graphic  d  splay  following  the  original  USAFSAM 
approach,  we  now  present  the  first  jvel  of  processing  by  our  matched  filter- 
approach. 

6.2 .2  The  Matched  Filter  Approach  or  the  Human  Subjects 

The  matched  filter  approach  requires  a  template  of  the  VER  and  a  mod'=l 
for  the  background  noise.  First,  w;  compute  a  raw  template  by  averaging 
over  the  data  sections  known  to  contain  the  VER.  Second,  the  fundamental 


reversal  rate  and  the  harmonics  (up  to  M5  Hz,  for  now)  are  picked  out  and 
used  to  obtain  the  smoothed  template.  For  the  background  noise,  we  rely  on 
the  previous  analysis  from  the  USAFSAM  epidural  recordings  (in  Macaca 
mulatta)  which  suggested  that  an  autoregressive  second-order  process  (al¬ 
though  of  variable  variance)  provided  a  good  approximation.  Through  a  few 
tost  runs,  we  found  some  good  values  for  the  present  scalp-electrode  record¬ 
ings.  These  fixed  values  were  used  throughout  our  analysis.  We  point  out, 
however,  that  apparently  no  such  model  is  needed  when  the  high-pass  filtered 
information  channels  are  used. 

The  results  of  using  the  matched  filter  for  the  same  data  sections  and 
the  center  electrode  (channel  #1) — as  just  described  for  the  FFT  method 
(except  that  precisely  40  sec  of  data  were  used) — are  shown  in  Figure  6.7, 
Plot  b,  for  the  female  subject  (labeled:  VFRDJA1),  and  in  Figure  6.8,  Plot  b, 
for  the  male  subject  (labeled:  VERKJ1 ) .  Observe  the  improvements,  espec¬ 
ially  for  the  male  subject.  As  for  the  FFT  method,  we  computed  standard 
deviation  SIG-Y  and  average  response  Y-3AF  for  the  individual  response  esti¬ 
mates  Y-HAT.  Numerical  values  are  given  in  Appendixes  2.1  -2.4. 

The  factors  in  the  figures  refe^  to  the  improvement  of  the  method  over 
the  FFT  method,  and  indicate  the  time  savings  which  can  be  achieved.  For 
example,  a  factor  of  10  implies  the  FFT  method  has  to  use  measurement  times 
10  times  as  long  to  achieve  the  same  accuracy  as  the  particular  alternative 
method.  Conversely,  by  our  method,  measurement  time  can  be  cut  to  one-tenth 
that  of  the  FFT  method  while  achieving  the  same  resolution.  A  summary  of 
the  standard  deviations  of  the  individual  estimates  is  shown  in  Table  6.1. 
From  these  values  the  improvement  factor  is  computed  by  taking  the  square  of 
the  ratio  of  corresponding  standard  deviations.  The  factors  in  the  figures 
include  the  5-percent  correction  factor,  which  is  due  to  having  used  5  per¬ 
cent  more  data  windows  in  the  FFT  method  than  in  the  corresponding  matched 
filter  analysis. 

Now  that  the  improvement  of  estimation  by  use  of  the  matched  filter  has 
been  seen,  the  next  interesting  step  is  to  combine  the  information  from  the 
individual  channels  by  suitable  regressior  . 

6.2.3  Integrated  Information  of  Response  Estimates 

Here,  a  suitable  approach  was  to  form  a  linear  combination  of  the  re¬ 
sponse  estimates  from  the  individual  matched  filters.  The  necessary  coeffi¬ 
cients  are  obtained  by  a  method  motivated  by  a  maximum  likelihood  argument 
and  regression,  respectively.  Although  (tor  our  current  understanding  of 
the  VER),  only  part  of  the  information  contained  in  the  outputs  of  the 
matched  filters  has  actually  been  used,  a  significant  further  reduction  of 
variability  could  be  achieved.  The  resul- s  from  combining  the  channels  1, 
2,  and  3  can  be  seen  In  the  displays:  in  Figure  6.7,  Plot  e,  for  the  female 
(labeled:  VERDJ123) ;  and  in  Figure  6.8,  pjot  c,  for  the  male  subject  ( la¬ 
beled:  VERKJ123).  Compare,  also,  the  standard  deviations  for  these  condi¬ 
tions  as  shown  in  Table  6.1.  Given  in  section  6.4  are  details  of  the 
estimation  of  the  weighting  factors  for  combining  outputs  from  the  matched 
filters.  Numerical  results  are  listed  In  Appendixes  3.1  -  3.4. 
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Figure  6.8:  Plots  ti  -  c.  Sequence  of  20  experimental  measurements 
for  a  male  subject .  Observe  the  variability  of  re¬ 
sponse  estimates. 


TABLE  6.1.  STANDARD  DEVIATIONS  OF  INDIVIDUAL  ESTIMATES 


Method 


No  stimulus 


Stimulus 


’ol lowing  this  rather  intuitive  explanation  of  our  approach,  we  now 
give  a  more  technical  presentation  and  evaluation  of  the  methods  we  devel¬ 
oped. 


6.3  TECHNICAL  ASPECTS  OF  MATCHED  FILTER  APPROACH  FOR  HUMAN  SUBJECT  DATA 


This  portion  of  the  report  deals  with  the  individual  programming 
steps  of  our  hierarchical  matched  filter  approach  (described  in  more  detail 
in  the  following  section).  Here,  we  give  a  brief  description  of  the  data 
processing,  with  some  examples  of  typical  results. 


We  begin  with  the  extraction  of  the  template,  turn  to  the  smoothing 
of  the  template,  review  the  matched  filter  and,  finally,  describe  the 
integration  of  the  information  from  tie  individual  data  channels.  A  summary 
of  the  individual  modeling  and  processing  stages  is  shown  in  Figure  6.6.  We 
also  include  a  description  of  the  resilts  shown  in  Figures  6.7  and  6.8,  and 
in  Table  6.1. 


6.3.1  Estimating  the  Template 


The  template  of  the  VER  is  obtained  in  a  two-step  procedure. 


First,  a  raw  estimate  is  obtained  by  averaging,  synchronously  with 
the  reference  channel,  a  1-sec  segment  of  data  over  all  data  which  are  known 
to  contain  the  VER;  e.g.,  for  each  individual,  the  VER  was  present  during 
about  T  =  60  sec  to  T  =  140  sec  in  each  of  the  10  repetitions  of  the  experi¬ 
ment.  Thus,  for  each  individual,  a  particular  raw  template  is  obtained.  These 
raw  templates  for  the  male  subject  arc.  displayed  in  Figure  6.9. 


Second,  this  template  is  smoothed.  The  smoothing  consists  of  a 
selection  of  the  fundamental  reversal  frequency  plus  some  of  its  harmonics. 
When  the  stimulus  frequency  is  15  Hz,  we  select  typically  only  the  15-Hz, 
30-Hz,  and  the  45-Hz  components  for  t  iis  smoothed  template.  Thus,  we  avoid 
difficulties  with  60-Hz  interference.  The  selection  is  obtained  by  finding 
the  FFT  of  the  raw  template,  zeroing  mt  the  undesired  frequencies,  and 
back-transforming  the  template  (Fig.  >.10).  Observe  the  asymmetry  of  these 
templates. 


We  have  not  tried  to  use  any  higaer  harmonics  at  this  time,  al¬ 
though  harmonics  up  to  about  100  Hz  may  be  present  in  VERs  (4).  Use  of 
higher  frequencies  may  possibly  contribute  to  further  reduction  of  the 
variability  of  VER  estimates. 
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6.3.2  The  Matched  Filter  for  the  H  iman  Subjects 


A  schematic  drawing  of  the  matched  filt'r  is  shown  in  Figure  6.11. 
Observe  that,  prior  to  the  matching,  the  raw  data  and  the  smoothed  template 
are  transformed  by  an  inverse  autor 3gres3ive  process  in  order  to  whiten 
the  background  of  the  EEG — the  noise  for  much  of  our  analysis.  The  need  for 
this  whitening  arises  from  the  fact  that  the  matched  filter  operates 
optimally  only  under  these  conditions. 

For  our  purposes,  we  found  that  autoregressive  model  values  ai  = 

1.6  and  03  =  ”.6  gave  good  results  regardless  of  the  particular  subject  or 
scalp  electrode.  These  values  resulted  in  a  slightly  stronger  high-pass 
filter  effect  than  those  we  had  derived  for  the  epidural  electrode 
recordings  (6).  Considering  the  low-pass  filter  characteristics  of  the 
transmission  of  VERs  from  the  visual  cortex  to  the  scalp,  these  new 
autoregressive  values  were  quite  reasonable,  even  though  they  were  found 
empirically.  We  also  noted  that,  for  the  data  which  were  recorded  by  means 
of  the  high-pass  filter  (Fig.  6.4),  no  preprocessing  with  the  inverse 
autoregressive  filter  was  needed  to  achieve  near-optimal  results.  This 
finding  was  anticipated  from  the  results  we  had  obtained  from  the  analysis 
of  earlier  recorded  data  (6).  For  m  example  of  a  typical  output  from  the 
matched  filter,  see  Figure  6.12. 

An  additional  benefit  may  be  obtained  by  use  of  the  matched  filter 
approach.  We  suspect  physiological  variability  in  the  amplitudes  of  the 
individual  evoked  frequencies  (e.g. ,  15  Hz  and  its  harmonics)  to  be 
omnipresent  but,  apparently,  either  in  an  independent  fashion  or  in  some 
counter-phase  fashion. 


6.4  INTEGRATION  OF  INFORMATION  FROM  INDIVIDUAL  INFORMATION  CHANNELS 

In  thi3  section,  our  objective  is  to  relate  the  estimated  VERs 
(from  the  individual  information  channels) to  the  stimulus  in  such  a  way  that 
the  stimulus  can  be  estimated  with  minimal  uncertainty.  We  may  regard  the 
condition  of  no  stimulus  to  correspond  to  y  «  0,  and,  when  the  30-percent 
contrast  stimulus  is  present,  to  y  -  1 ,  the  "unit  stimulus."  All  other 
potential  stimuli  would  be  scaled,  in  an  appropriate  way,  relative  to  these 
two  experimental  conditions. 

Given  VER  estimates,  x, we  nay  then  ask  how  to  combine  them  to 
estimate  a  stimulus  condition  (or  some  "generalized  response")  —  e.g.,  y — 
with  least  variability.  More  forma  . Ly,  we  wish  to  find  a  method  to  obtain. 

min  [d(y-y)/xi ,  xg,...,  Xfji,  (6.1) 

where  d  is  some  suitable  metric.  Vc  are  mainly  interested  in  the  usual 
least  squares  sense;  that  is, 

d(y-y)  -  (y-y)2,  (6.?) 


SC 


RMP 


I .1 


which  looks  very  much  like  the  usual  regression  problem.  However,  here  one 
has  to  be  careful  in  obtaining  weights  Bi,...,6n;  for  the  errors  in  the 
model  do  not  lie  in  the  stimulus  conditior,  y  (which  is  known  from  the 
experiment),  as  in  the  usual  regression  model: 


y  =  Bo  +  Mi  + 


Bnxn  +  e 


(6.3) 


but  in  the  x's,  as  in: 


xi  =  Yiy  +  ei 


i=1 _ n  . 


(6.M) 


This  situation  is  known  as  the  "error  invariable"  problem.  In  many 
situations,  the  usual  regression  estimation  for  the  parameters  Yj,...,Yn  is 
still  useful.  In  fact,  it  also  leads  to  good  estimation  for  our  particular 
problem.  Thus,  we  may  obtain  regression  coefficients  by  use  of: 


y  =  [x x  x  r1  x' 


(6.5) 


where,  as  usual,  the  rows  of  X  are  formed  by  "x"  for  the  individual 
experiments. 


Still,  a  more  rigorous  argument  about  the  use  of  regression 
coefficients  (or,  equivalently,  weight  coefficients)  is  desirable.  To 
clarify  the  situation,  refer  to  the  displays  in  Figure  6.1 3.  The  regression 
—  which  is  fit  to  such  data,  minimizing  ’•.he  residual  error  in  y  —  will 
tend  to  be  too  flat.  Thus,  when  y  is  predicted  on  the  basis  of  observations 
x,  using  the  regression  line,  biased  estimates  will  result. 

In  this  particular  case,  two  remedies  are  apparent:  one  appears  to 
bn  a  quick  fix;  the  other  is  motivated  by  maximum  likelihood  estimation.  As 
can  easily  be  shown,  they  are  identical  here. 


The  quick  fix,  which  seems  to  be  to  pull  the  regression  line 
through  the  origin  of  the  y  and  x  space,  nay  be  done  either  by  fitting  a 
model  without  the  regression  constant  60  or,  alternatively,  by  adding  a  lot 
or  pseudo  data  points  at  the  origin  while  fitting  a  model.  This  latter 
approach  is  convenient  when  standard  regression  packages,  which  often 
include  80 »  are  to  be  used. 
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How  error  in  the  assumed  independent  variable,  x,  yields 
a  biased  regression  line  for  assumed  dependent  variable, 
(which  will  usually  be  identified  with  the  stimulus 
trength) . 


With  regard  to  the  maximum  likelihood  (ML)  estimation  (MLE) , 
consider  the  following  simple  model  for  n= correlated  channels  (Fig.  6.1H) 
for  two  stimulus  conditions  y  =  0  and  y  =  1  ; 


*i  =  "Y i  Yj  +  Gi;  i=1  •••  n,  j-0,1,  (6.6) 


and 


COV  [£,£]»  C,  (6.7) 

where  Yi  stands  for  the  fixed  part  of  the  relation  between  xj  and  y,  while 
Gj  stands  for  the  random  part  of  that  relationship.  Then, 

p ( x  /  y,  Y,  C)  =  c  exp  {-1/2  (x  -  Y  y)/  C  -1(x  -  Y  y).  (6.8) 


The  maximum  likelihood  estimate  of  Y.>  given  the  stimulus  present  condition 
(y  =  1 ) ,  is  simply: 

Y  =  x  /  y  ,  (6.9) 


where  the  averaging  is  performed  over  independent  (for  us,  successive) 
experiments.  When  these  gain  factors,  summarized  in  _Y,  have  been  obtained, 
one  may  conversely  proceed  to  estimate  y,  given  x.  The  maximum  likelihood 
estimate  for  y  is  by  equation  6.8,  and  by  replacing  the  estimate  _Y  for  y_: 

y  ~  y  C-1x  /  Y  C'1  Y  .  (6.10) 


We  can  show,  fairly  simply,  that  estimation,  via  equation  6.10 
leads  to  the  same  estimates  as  use  of  the  regression  formula,  equation  6.5, 
when  an  infinitely  large  number  of  points  at  the  origin  is  added.  This 
finding  is  due  to  the  equivalence  (up  to  a  multiplicative  factor,  in  this 
case)  of  x'^withj,  and  (x'x)  with  C. 

We  should  point  out  at  this  time  that  the  model  of  the  estimates  of 
VERs,  as  shown  in  equation  6.7,  is  not  quite  realistic  in  view  of  our 
current  data  analysis:  that  C  changes,  generally  speaking,  increases 
monotonic  with  stimulus  intensity  (Fig.  6.15).  We  presume,  also,  that  the 
contours  of  the  distribution  of  estimates  change  with  stimulus  intensity. 
Some  preliminary  investigations  along  these  lines  showed  the  usefulness  of  a 
variability  measure  for  predicting  stimulus  condition.  Since  these  aspects 
exceed  by  far  the  scope  of  our  current  effort,  however,  we  have  to  leave 
further  investigation  for  the  future. 


In  summary,  significant  reduction  of  the  estimation  of  VERs  is 
possible  by  a  combined  matched  filter  and  maximum  likelihood  approach. 
However,  the  current  approach  is  hierarchical  and  should  not  be  presumed  to 
be  optimal  with  regard  to  potential  reduction  of  variability.  One  of  the 
important  questions  is  how  one  might  go  about  finding  an  integrated  (not 
hierarchical)  approach  which  has  the  potential  to  yield  near  optimal 
performance.  The  necessary  approach  seems  to  be  dictated  by  our  improved 
understanding  of  the  VER.  Some  of  the  important  findings,  which  will  lead 
eventually  to  near  optimal  assessment  of  visual  performance,  are  discussed 
in  the  following  pages. 
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r.  SOME  BASIC  STUDIES  AND  PROPERTIES  OF  THE  HUMAN  VER 


Three  properties  of  the  data  appear  particularly  interesting  for 
refining  our  understanding  of  the  VER.  They  concern:  the  variability  of 
the  VER;  the  mutual  dependencies  of  harmonic  frequencies  and  fundamental 
frequency;  and  the  effect  of  blinking  on  the  phase  angle  of  the  VER.  The 
observations  are  in  accordance  with  our  understanding  of  the  underlying 
physiology  of  information  processing,  and  are  discussed  here  in  this 
context. 


7.1  THE  CONCEPT  OF  VARIABILITY 

As  already  mentioned  (Figs.  6.13  and  6.15),  the  variability  of 
response  estimates  in  individual  channels  (corresponding  to  differential 
voltages  between  electrodes) increases  when  a  stimulus  is  presented.  This 
finding  is  true  for  both  individual  channels.  The  increase  is  strong  enough 
to  be  useful  for  predicting  stimulus  condition. 

In  general,  we  computed,  for  individual  40-sec  segments  of  data, 
the  standard  deviation  of  individual  VER  estimates  based  on  10-sec 
subsegments.  Thus,  the  variability  of  4  subsegments  generated  a  variability 
index.  This  index  could  be  used,  in  turn,  in  the  (already  described)  maximum 
likelihood  estimation  for  assessing  stimulus  condition.  In  fact,  we  found 
from  regression  analysis  that  some  50  -  70  percent  of  the  regression  model 
was  explained  by  this  variable. 

Conceptually,  what  appears  to  happen  is  that  the  "gain"  of  the 
visual  system  varies  in  a  random  fashion.  Thus,  the  larger  an  input  signal 
(e.g. ,  no-contrast  »  0  vs.  unit-contrast  =1),  the  larger  will  be  the 
absolute  variation  of  the  response.  Thus,  such  variation  relates  to  the 
amplitude  of  the  stimulus  intensity. 

A  second  aspect  was  highly  interesting  to  us.  It  concerns  the 
hemispheric  differences  of  the  VER.  We  found  that,  for  the  female  subject, 
mainly  the  left  hemisphere  gave  large  variability  —  while  the  right 
hemisphere  delivered  a  much  more  constant  response.  In  the  male  subject, 
the  converse  was  true:  the  right  hemisphere  gave  highly  variable  responses 
—  while  the  left  hemisphere  provided  a  very  constant  response.  This  result; 
which  may  be  a  pure  fluke  of  the  "random  sampling  of  the  population"  (not  of 
our  data  processing),  is  in  accordance  with  the  general  outline  in  our  final 
report  for  phase  1(5)  .There,  we  discussed  the  hemispheric  differences,  in 
accord  with  Beaumont  et  al.  (1),  Ben-Dov  and  Carmon  (2),  and  others,  the 
female  utilizes  her  left  hemisphere  more, as  compared  with  the  male;  and, 
conversely,  the  male  utilizes  more  of  his  right  hemisphere  than  does  the 
female.  Thus,  in  the  male,  die  right  hemisphere  appears  to  be  more  complex; 
and,  in  the  female,  the  left  hemisphere,  associated  with  the  differing 
functions  of  these  hemispheres.  This  finding  is  in  accordance  with  our 
observation  and  has  an  obvious  bearing  on  what  might  be  expected  from  the 
VER  in  these  hemispheres — the  VER  m  the  more  complex  hemisphere  is  expected 
to  be  more  variable.  Substantiating  this  finding  with  more  representati ve 
random  samples  of  the  general  populat  Lon  will  be  important . 


The  increased  variability  was  reflected  in  a  much  smaller  weight  factor 
for  that  hemisphere  when  the  individual  channels  were  combined  in  the  maxi¬ 
mum  likelihood  estimation.  In  fact,  a  slightly  negative  weight  factor  was 
estimated,  both  for  the  female  and  the  male.  The  fact  that  the  coefficients 
are  not  significantly  negative  might  be  an  indication  of  a  counter-phase 
variation  of  the  VER  between  the  hemispheres:  when  one  hemisphere  receives 
more  input,  the  other  one  receives  less.  Such  a  model  has  to  be  considered 
in  view  of  the  known  existence  of  competing  neural  networks  (mutual  inhibi¬ 
tory  innervation). 

7.2.  RELATION  OF  FUNDAMENTAL  AND  HARMONIC  FREQUENCIES 


This  study  was  of  limited  nature,  and  was  performed  only  in  the  female 
because  of  her  strong  and  clear  VERs.  Basically,  we  tried  to  determine  any 
systematic  relation  in  the  variation  of  the  fundamental  and  harmonic  variabil¬ 
ity  of  response  estimates  (or  better  estimators  of  stimulus  condition)  could 
be  reduced.  For  this  purpose,  we  checked  the  autocorrelation  of  15-Hz  ampli¬ 
tude  estimates  from  successive  1-sec  data  windows  with  the  corresponding 
amplitude  estimates  at  30  Hz  (Fig.  7.1.).  From  a  statistical  viewpoint,  no 
significant  correlation  could  be  observed,  for  example,  on  a  95-percent  level. 
Thus,  for  the  present,  we  disregarded  any  potential  relation. 

In  retrospect,  the  situation  may  be  different  for  the  male  subject, 
because  variability  of  the  "true"  VER  seems  to  be  a  much  more  significant 
component  when  compared  with  the  variability  induced  by  the  noisy  background 
EEC.  Thus,  a  possible  relationship  may  become  statistically  more  signifi¬ 
cant. 

7.3  THE  EFFECT  OF  BLINKING 

For  some  time,  blinking  has  been  our  concern  as  one  of  the  contributing 
factors  in  variability  of  response  estimates  in  human  subjects.  Thus,  we 
performed  an  additional  experiment,  in  which  the  individual  was  asked  to 
blink,  self-controlled,  every  5  sec.  Accuracy  was  verified  from  a  separate 
EMC  signal  channel  (described  earlier) . 

The  major  finding  of  this  limited  effort  was  the  change  of  the  phase 
angle  of  the  VER  during  blinking.  For  both  subjects,  the  phase  angle  was 
reduced  (for  the  fundamental  reversal  rate  and,  in  corresponding  fashion, 
for  the  harmonics).  This  reduction  in  phase  angle  correspond-  to  a  faster 
response  to  the  stimulus.  For  the  male  subject,  the  change  corresponded  to 
some  7  ms;  and,  in  the  female,  about  1.2  ms.  Although  these  changes  may 
appear  small,  they  are  significant. 

We  came  to  this  conclusion  in  the  following  way:  We  computed,  in 
synchrony  with  the  reference  channel,  the  FFT  of  successive  1-sec  complex 
planes.  In  the  event  that  no  stimulus  was  present,  a  random  walk  of  this 
sum  developed  (Fig.  7.2).  On  the  other  hand,  when  a  stimulus  is  present,  a 
consistent  component  will  be  in  one  particular  direction.  These  results  are 
shown  for  both  subjects  in  Figures  7.3  and  7.4. 


59 


riots  a  ana  _D.  comparison  of  phase  angle  of  VER  with  and  without  blinking. 
Blinking  reduces  the  phase  angle.  The  corresponding  change  of  delay  time  is 
1.2  ms.  Subject:  female;  15-Hz  reversal  rate. 


The  left-hand  side  shows  the  results  of  integrating  the  complex  vectors 
when  the  subject  suppressed  blinking  (reduced  to  some  2  to  3  times/min), 
while  the  right-hand  side  displays  the  result  for  regular  blinking  every  5 
sec.  For  both  subjects,  a  drastic  reduction  of  the  phase  angle  results. 
(When  interpreting  the  figures,  observe  the  unequal  scales  in  the  left-hand 
part.)  The  result  is  interesting,  especially  because  one  will  usually  asso¬ 
ciate  late  responses  with  fatigue  and,possibly*more  variability. 

In  any  event,  for  reducing  variability  in  response  estimates  in  the 
presence  of  blinking,  the  incorporation  of  the  effect  of  blinking  on  the 
phase  lag  of  the  VER  seems  to  be  important.  The  importance  is  more  striking 
when  one  considers  higher  harmonics.  There,  changes  of  phase  angle  in 
excess  of  90  degrees  are  observed.  Clearly,  from  the  properties  of  the 
matched  filter,  this  knowledge  about  phase  shifts  is  crucial. 

At  this  time,  however,  we  could  not  find  any  indication  for  reducing 
variability  of  response  estimates  by  incorporating  phase  variation  into 
estimation  of  the  VER  when  the  subject  is  not  blinking.  Over  periods  of 
more  than  2  min,  in  the  present  recordings,  blinking  seemed  to  occur  rarely 
(e.g.,  3  or  4  times).  Due  to  our  doubt  about  the  accuracy  of  our  channel 
for  recording  blinking,  we  decided  to  verify  performance  of  that  channel  by 
encouraging  the  subject  to  blink  regularly  every  5  sec.  We  can  only  suspect 
that,  when  the  subject  is  encouraged  to  blink  freely,  the  consideration  of 
phase  lags  becomes  important  for  modeling  the  VER.  Inducing  a  subject  to 
blink  freely  may  be  attractive  because  of  the  aforementioned  faster  re¬ 
sponse,  which  may  be  an  expression  of  increased  alertness. 

In  summary,  we  made  a  number  of  interesting  observations  with  regard  to 
the  variability  of  the  VER  and  its  causes.  Because  of  time  constraints, 
however,  a  number  of  questions  could  not  be  thoroughly  investigated.  These 
had  to  be  left  to  future  analysis  and  to  new  experimental  setups. 


7.^  STUDY  OF  SEGMENTING 

Following  our  observation  of  the  possible  clusters  of  the  variance  of 
the  EEG  background  activity  in  data  of  Macaca  mulatta  (Fig.  4.12),  we  tested 
the  significance  of  weighing  data  according  to  the  residual  power  in  the 
matched  filter  (the  sum  of  squares  of  the  residuals  within  1-sec  segments 
of  data).  The  following  two  opposing  models  were  tried — 

First:  The  relative  weight  of  the  contribution  of  an  amplitude 
estimate  was  reduced  when  the  residual  power  increased.  This  processing 
counteracts  a  model  in  which  corrupt, ing  noise  is  assumed  to  be  superimposed 
on  the  VER.  Such  a  processing  method  is  attractive,  because  it  leads 
naturally  to  a  robust  procedure  (one  which  is  not  too  badly  affected  by 
sudden  artifacts,  such  as  possible  loosening  of  electrodes).  Second:  We 
tested  a  processing  which  increases  the  weight  of  amplitude  estimates  within 
these  1 -sec  data  windows  when  the  residual  power  increases.  Interestingly, 


I 
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we  did  not  find  any  statistically  significant  differences  in  these  two 
alternative  processing  methods  compared  to  the  plain  integration  in  time  of 
the  matched  filter  output.  We  assume,  therefore,  that  the  current  method 
must  be  in  a  sense  close  to  optimal.  However,  one  should  keep  in  mind  that 
this  attempt  was  not  based  on  any  of  the  (in  this  respect)  possibly  varying 
properties  of  the  human  subject  data.  Following  a  more  rigorous  analysis,  a 
more  sophisticated  data  processing  might  improve  estimation.  This  question 
has,  however,  to  be  left  to  future  investigation. 
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DISCUSSION 


According  to  the  results  we  have  obtained  to  date,  our  approach  appears 
to  be  a  valid  way  to  reduce  the  variability  of  response  estimates.  Our  anal¬ 
ysis  had  impact:  on  the  selection  of  the  differentials  between  electrodes 
which  were  to  be  used  as  information  channels;  on  the  prewhitening  of  record¬ 
ed  data  by  a  high-pass  filtering  process  (eliminating  the  problem  of  ampli¬ 
fier  saturation,  and  simplifying  subsequently  the  matched  filter  which  does 
not  require  the  inverse  autoregressive  filter);  and  on  the  ultimate  structure 
of  our  hierarchical  processing  method. 

The  present  data  indicate  that  hemispheric  differences — possibly  linked 
to  the  sex  of  the  subject — and  the  effect  of  blinking  may  be  some  of  the 
next  variables  to  be  investigated.  Other  variables  are  day-to-day 
variations  which  may,  in  part,  be  due  to  variation  in  electrode  locations. 

A  word  should  also  be  said  about  the  current  data-acquisition  method. 

We  found  the  time  resolution  to  be  somewhat  of  a  limiting  factor  because  of: 

(a)  considerable  random  tape- speed  variations  up  to  10  ms/sec;  and 

(b)  a  sampling  rate  of  only  256  samples/sec. 

This  low  sampling  rate  results,  in  view  of  the  random  tape-speed  variations, 
in  a  large  uncertainty  in  phase  angles  at  higher  harmonics  of  the  VER. 

The  problem  can  be  cured,  to  some  extent,  by  more  sophisticated 
software  which  synchronizes  the  data  processing  with  the  reference  channel 
more  often  than  once  per  second.  Possibly  a  simple  model  of  the 
tape-speed-sensitive  FM  decoder  of  the  tape  deck  should  be  included,  and 
interpolation  of  data  and  reference  signal  between  sample  points  (for  proper 
phase  alignment  relative  to  previous  data  sections)  should  be  performed.  Of 
course,  if  data  acquisition  did  not  have  to  go  through  this  replay  procedure 
with  the  tape  machine  (as  with  the  use  of  a  separate  computer  for  data 
acquisition),  these  problems  would  disappear  automatically.  Hence,  the 
software  would  remain  much  simpler. 

Beyond  the  current  methodology,  we  feel  that  importance  lies  in 
developing  more  rigorous  dynamic  models  of  the  VER  (not  just  a  fixed  matched 
filter) ,  which  make  explicit  use  of  the  variability  of  the  VER.  The  use  of 
different  stimulus  conditions — possibly  a  superposition  of  two  reversal 
rates  at,  e.g.,  14  Hz  and  17  Hz — may  permit  the  identification  of  mechanisms 
which  contribute  to  the  variability  of  the  VER.  When  a  model  for 
variability  of  components  and  their  mutual  relation  can  be  developed, 
further  reduction  of  uncertainty  associated  with  the  VER  can  be  expected. 

From  a  philosophical  point  of  view,  feeding  a  more  complex  stimulus  into  the 
visual  system  (viz.,  14  Hz  and  17  Hz)  might  permit  us  to  extract  more 
information  from  the  VER  (provided  the  channel  capacity  of  the  system  is  not 
exceeded) ,  and  to  surpass  the  results  we  were  able  to  achieve  so  far  in 
terms  of  reduced  variability  and/or  necessary  experimental  measurement  time. 
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APPENDIX  1.1:  STATISTICS  FOR  FFT  METHOD,  CENTER  ELECTRODE 
NO  STIMULUS— FOR  FEMALE  SUBJECT 
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STATISTICS  FOR  FFT  METHOD,  CENTER  ELECTRODE, 
WITH  STIMULUS— FOR  FEMALE  SUBJECT 
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APPENDIX  1.3:  STATISTICS  FOR  FFT  METHOD,  CENTER  ELECTRODE 
NO  STIMULUS — FOR  MALE  SUBJECT 
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0. 1791E+01  0. 1000E+00 

0. 1549E+01  0 . 1000E+01 

0.2016E+01  0 . 1000E+01 

0.2309E+01  0.1000E+01 

0 . 1937E+01  0. 1000E+01 

0.2511E+O1  0 . 1000E+01 

0 . 2592E+01  0.1000E+01 

0.2167E+01  0. 1000E+01 

0.2611E+01  0.1000E+01 

0. 1477E+01  0 . 1000E+0 1 

GIVE  STIMULUS  CONDITION,  E.G.  Y  -  1, 

0 

REG-CONSTANT  +  COEFFICIENTS 
-0.14385505  0.54573237 


Y-HAT 

* 

-0.1439 

+ 

0.5457  XI  + 

N'SAMP 

- 

1 

Y-HAT 

0.0213 

Y-BAR 

-  0.0213 SIG-Y 

-  0.0000 

Y-SHARP 

0.0000 

NSAMP 

2 

Y-HAT 

-0.0965 

Y-BAR 

-  -0.376  SIG-Y 

-  0.0832 

Y-SHARP 

-2.2141 

NSAMP 

* 

3 

Y-HAT 

0.1441 

Y-BAR 

-  0.0230  SIG-Y 

-  0.1203 

Y-SHARP 

=» 

5.2365 

NSAMP 

» 

4 

Y-HAT 

0.0223 

Y-BAR 

-  0.0228 SIG-Y 

-  0.0982 

Y-SHARP 

4.3072 

NSAMP 

5 

Y-HAT 

0.0250 

Y-BAR 

-  0.233  SIG-Y 

-  0.0851 

Y-SHARP 

a 

3.6584 

NSAMP 

» 

6 

Y-HAT 

-0.0374 

Y-BAR 

-  0.0131  SIG-Y 

-  0.0800 

Y-SHARP 

* 

6.0894 

nsamp 

- 

7 

Y-HAT 

-0.0146 

Y-BAR 

-  0.0092 SIG-Y 

-  0.0738 

Y-SHARP 

Si 

8.0446 

NSAMP 

- 

8 

Y-HAT 

-0.0903 

Y-BAR 

-  -0.0033 SIC-Y 

-  0.0768 

Y-SHARP 

a 

-23.5287 

NSAMP 

- 

9 

Y-HAT 

0.0386 

Y-BAR 

-  0.0014  SIG-Y 

-  0.0732 

Y-SHARP 

a 

52.7674 

nsamp 

- 

10 

Y-HAT 

-0.0126 

Y-BAR 

-  0.0000  SIG-Y 

-  0.0692 

Y-SHARP 

a 

-4768.6792 

A 

y 

current 
estimates 
plotted  in 
Figure  6.8 
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APPENDIX  1.4:  STATISTICS  TOR  FFT  METHOD,  CENTER  ELECTRODE, 
WITH  STIMULUS — FOR  MALE  SUBJECT 


$  R  REGHOD 
La  or  DATA  FILE 
90 

#  OP  INDEPEDENT  VARIABLES ,E.G.  3 

1 

#  OF  DIFFERENT  STIMULUS  CONDITIONS ,E.G.  2 

2 

#  OF  POINTS  IN  GROUP,  E.G.  10 

10 

KJ  t 

KJ  1 


x  y 


0.3026E+00 

O.OOOOE+OO 

0.8686E-01 

0.0000E+00 

0.5277E+00 

O.OOOOE+OO 

0.3045E+00 

O.OOOOE+OO 

0.3095E+00 

O.OOOOE+OO 

0. 1950E+00 

O.OOOOE+OO 

0.2368E+00 

O.OOOOE+OO 

0.9807E-01 

O.OOOOE+OO 

0.3344E+00 

O.OOOOE+OO 

0.2404E+00 

O.OOOOE+OO 

0. 1791E+01 

0. 1000E+00 

0. 1549E+01 

0. 1000E+01 

0 . 2016E+01 

0 .  lOOOE+Ol 

0 . 2309E+01 

0 . 1000E+0 1 

0. 1937E+01 

0. 1000E+01 

0.251 1E+01 

O.IOOOE+Ol 

0.2592E+01 

0. 1000E+01 

0.2167E+01 

0. lOQOE+Ol 

0.761 1E+01 

0 . 1000E+01 

0 . 1477E+0 1 

0. lOOOE+Ol 

STIMULUS  CONDITION.  E.G.  Y  - 

1 

REG-CONSTANT  +  COEFFICIENTS 
-0.14385505  0.54573237 


selected 


Y-HAT  =  -0,1439'+  0.5457  XI  + 


NSAMP 

* 

1 

Y -1IAT 

* 

0.8337 

Y-BAK 

= 

0.8337 

SIG-Y 

-  0.0000 

Y-SHARP 

- 

0.0000 

NS  AMP 

= 

2 

y-HAi 

* 

0.7017 

Y-BAK 

* 

0.7677 

SIG-Y 

*  0.0934 

Y -SHARP 

0.1216 

NSAMP 

= 

3 

Y-HAT 

= 

0.9565 

Y-BAK 

0.8306 

SIC-Y 

*  0.1274 

Y-SHARP 

* 

0.1534 

to  AMP 

4 

Y-HAT 

1.1161 

Y-BAK 

0.9020 

SIG-Y 

*  1.766 

Y-SHARP 

- 

0.1958 

NSAMP 

= 

3 

Y-HAT 

0.9130 

Y-BAK 

0.9042 

SIG-Y 

-  0.1530 

Y-SHARP 

- 

0.1693 

NS  AMP 

» 

f) 

Y-HAT 

1.2268 

Y-BAK 

0.9580 

SIG-Y 

-  0.1899 

Y-SHARP 

* 

0.1983 

NSAMP 

* 

7 

Y-HAT 

* 

1.2708 

Y-BAK 

1.0027 

SIC-Y 

-  0.2099 

Y-SHARP 

- 

0.2093 

NSAMP 

- 

8 

Y-HAT 

m 

1.0386 

Y-BAK 

1.0072 

SIC-Y 

-  0.1947 

Y-SHARP 

* 

0.1933 

NSAMP 

• 

9 

Y-HAT 

m 

1.2811 

Y-BAK 

1.0376 

SIG-Y 

»  0.2037 

Y-SHARP 

m 

0.1964 

NSAMP 

- 

10 

Y-HAT 

“ 

0.6619 

Y-BAK 

1.0000 

SIO-Y 

*  0.2259 

Y-SHARP 

- 

0.2259 

A 

y 

current 
estimates 
plotted  in 
Figure  6.8 
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L"  HI  -  LF 


APPENDIX  2.1:  STATISTICS  FOR  MATCHED  FILTER  METHOD,  CENTER 
ELECTRODE,  NO  STIMULUS — FOR  FEMALE  SUBJECT 


LUROREGMOD  FILE 
90 

II  OF  INDEPENDENT  VARIABLES, E.G.  3 
3 

II  OF  DIFFERENT  STIMULUS  CONDITIONS,  E.G. 
2 

II  OF  POINTS  IN  GROUP, E.G.  10 
10 

TEMP ID. OUT  1 

TEMP2D.0UT  2 

TEMP3D.0UT  3 

TEMP ID. OUT  1 

TEMP2D . OUT  2 

TEMP3D.0UT  3 


VERDJ1 , 


stimulus  condition 


X1 

x2 

X3 

y 

0. 9931E-03 

0. 6960E-01 

0. 3915E-01 

O.OOOOE+OO 

0. 3194E-01 

0. 6296E-01 

0. 1829E-02 

0 . OOOOE+OO 

0. 6406E-01 

-0.8106E-02 

-0.1875E+00 

O.OOOOE+OO 

0. 2677E-02 

0.2856E-01 

0. 2438E+00 

O.OOOOE+OO 

0. 3547E-01 

0. 1780E-01 

-0.1254E+00 

O.OOOOE+OO 

0. 3902E-01 

-0.3840E-01 

-0. 4314E-01 

O.OOOOE+OO 

0. 1096E+00 

0. 6040E-01 

-0.8432E-01 

O.OOOOE+OO 

0. 4794E-01 

0.3290E-01 

-0. 6859E-01 

O.OOOOE+OO 

0.467 7 E-01 

0. 1066E+00 

0.1987E+00 

O.OOOOE+OO 

0.1385E-01 

0. 5961E-01 

0. 9772E-01 

O.OOOOE+OO 

0.1941E+01 

0. 2001E+01 

0. 1639E+01 

0. 1000E+01 

0. 2072E+01 

0. 2103E+01 

0. 1995E+01 

O.IOOOE+Ol 

0.1939E+01 

0.2029E+01 

0.1860E+01 

0. 1000E+01 

0.2014E+01 

0.2056E+01 

0.1869E+01 

O.IOOOE+Ol 

0.2068E+01 

0.2039E+01 

0.1986E+01 

0.1000E+01 

0. 1826E+01 

0. 1993E+01 

0. 1806E+01 

0. 1000E+01 

0. 2012E+01 

0. 1916E+01 

0.1671E+01 

0, 1000E+01 

0.1987E+01 

0. 1933E+01 

0. 2117E+01 

0. 1000E+01 

0.1934E+01 

0.1965E+01 

0.2377E+01 

0. 1000E+01 

0.2024E+01 

0. 2081E+01 

0. 2496E+01 

0. 1000E+01 

selected 


the  x  's  are  the 
endpoint  estimates  from 
the  matched  filter  output 


GIVE  STIMULUS  CONDITION,  E.G. 

0 

REG-CONSTANT  +  COEFFICIENTS 
-.02019  .5148  0  0 


Y-l 


negl ected 


Y-HAT 

-0.0202 

+ 

0.5148  XI  + 

0.0000  X2 

+  0. 

0000  X3  + 

N3AMP 

1 

Y-HAT 

-0.0197 

Y-8AR 

-0.0197 

S1G-Y 

-0.0000 

Y-SHARP 

0.0000 

NSAMP 

2 

Y-HAT 

-0.0037 

Y-BAR 

-0.0117 

SIC.-Y 

-0.0113 

Y-SHARP 

-0.9618 

NS, VIP 

=» 

3 

Y-HAT 

0.0128 

Y-BAR 

-0,0035 

S1G-Y 

•0.0162 

Y-SHARP 

-4.5775 

NSAMP 

4 

Y-HAT 

-0.0188 

Y-BAR 

-0.0074 

SIC-Y 

-0.0153 

Y-SHARP 

-2.0774 

NSAMP 

5 

Y-HAT 

-0.0019 

Y-BAR 

-0.0063 

SIG-Y 

-0.0135 

Y-SHARP 

-2.1459 

NSAMP 

6 

Y-HAT 

-0.0001 

Y-BAR 

-0.0052 

SIG-Y 

-0.0123 

Y-SHARP 

-2.3454 

NSAMP 

7 

Y-HAT 

0.0362 

Y-BAR 

0.0007 

SIG-Y 

-0.0193 

Y-SHARP 

28.3809 

NSAMP 

8 

Y-HAT 

0.0045 

Y-BAR 

0.0012 

SIC-Y 

-0.0179 

Y-SHARP 

15.4908 

NSAMP 

9 

Y-HAT 

0.0039 

Y-BAR 

0.0015 

SIC-Y 

*0.0168 

Y-SHARP 

11.4932 

NSAMP 

10 

Y-HAT 

-0.0131 

Y-BAR 

0.0000 

SIG-Y 

-0.0165 

Y-SHARP 

2069.2251 

y  given 
plotted 


in  Figure  6.7 
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APPENDIX  2.2:  STATISTICS  FOR  MATCHED  FILTER  METHOD,  CENTER 
ELECTRODE,  WITH  STIMULUS — FOR  FEMALE  SUBJECT 


R  REGMOD 
LU  OF  DATA  FILE 
90 

#  OF  INDEPENDENT  VARIABLES,  E.G.  3 
3 

#  OF  DIFFERENT  STIMULUS  CONDITIONS,  E.G. 


»  OF  POINTS  IN  CROUP,  E.C.  10 
10 

TEMP1D.0UT 
TEMP2D.0UT 
TEMP3D.0UT 
TEMP1D.0UT 
TEMP2D . OUT 
TEMP3D.0UT 


X1 

x2 

x3 

y 

0. 9931E-03 

0. 6960E-01 

0. 3915E-01 

) . OOOOE+OO 

0. 3194E-01 

0. 6296E-01 

0.1829E-02 

0 . OOOOE+OO 

0. 6406E-01 

-0.8106E-02 

-0. 1875E+00 

0. 0000E+00 

0.2677E-02 

0.2856E-01 

0. 2438E+00 

o.ooooe+oo 

0. 3547E-01 

0.1780E-01 

-0.1254E+00 

0. OOOOE+OO 

0.3902E-01 

-0.3840E-01 

-0.4314E-01 

0. OOOOE+OO 

0.1096E+00 

0.6040E-01 

-0.8432E-01 

O.OOOOE+OO 

0.4794E-01 

0.3290E-01 

-0. 6859E-01 

O.OOOOE+OO 

0.4677E-01 

0. 1066E+00 

0.1987E+00 

U.OOOOC+OO 

0. 1385E-01 

0. 5961E-01 

0. 9772E-01 

O.OOOOE+OO 

0. 1941E+01 

0. 2001E+01 

0. 1639E+01 

0. lOOOE+Ol 

0.2072E+01 

0.2103E+01 

0.1995E+01 

0. 1000E+01 

0. 1939E+01 

0. 2029E+01 

0. 1860E+01 

0 . 1000E+01 

0.2014E+01 

0.2056E+01 

0. 1869E+01 

0. 1000E+01 

0.2068E+01 

0.2039E+01 

0. 1986E+01 

0. lOOOE+Ol 

0. 1826E+01 

0. 1993E+01 

0. 1806E+01 

0. 1000E+01 

0. 2012E+01 

0. 1916E+01 

0.1671E+01 

0.1 OOOE+Ol 

0.1987E+01 

0. 1933E+01 

0. 2117E+01 

0. lOOOE+Ol 

0.1934E+01 

0.1965E+01 

0. 2377E+01 

0. lOOOE+Ol 

0.2024E+01 

0.2081E+01 

0.2496E+01 

0. 1 OOOE+Ol 

GIVE  STIMULUS  CONDITION, 

,  E.G.  Y-l 

the  x  'a  are  the 
endpoint  estimates  from 
the  matched  filter  output 


REG-CONSTANT  +  COEFFICIENTS 
-.02019  .5148  0  0 


neglected 


Y-HAT 

-0.0202  +  0.5148  XI  + 

'0.0000  X2 

+  0. 

0000  X3  + 

NSAMP 

1 

Y-HAT  -0.9791 

Y-BAR 

*0.9791 

SIG-Y 

-0.0000 

Y-SHARP 

0.0000 

NSAMP 

2 

Y-HAT  -1.0465 

Y-BAR 

-1.0128 

SIC-Y 

-0.0476 

Y-SHARP 

0.0470 

NSAMP 

3 

Y-HAT  -0.9778 

Y-BAR 

-1.0012 

SIG-Y 

-0.0393 

Y-SHARP 

0.0392 

NSAMP 

4 

Y-HAT  -1.0166 

Y-BAR 

-1 . 0050 

SIG-Y 

-0.0330 

Y-SHARP 

0.0328 

NSAMP 

5 

Y-HAT  -1.0445 

Y-BAR 

-1.0129 

SIG-Y 

-0.0336 

Y-SHARP 

0.0332 

NSAMP 

6 

Y-HAT  -0.9199 

Y-BAR 

*0.9974 

SIG-Y 

-0.0484 

Y-SHARP 

w 

0.0486 

NSAMP 

7 

Y-HAT  -1.0154 

Y-BAR 

-1 . 0000 

SIG-Y 

-0.0447 

Y-SHARP 

a* 

0.0447 

NSAMP 

8 

Y-HAT  -1.0028 

Y-BAR 

-1.0003 

SIG-Y 

-0.0414 

Y-SHARP 

m 

0.0414 

NSAMP 

9 

Y-UAT  -0.9756 

Y-BAR 

*0.9976 

SIG-Y 

-0.0396 

Y-SHARP 

m 

0.0397 

NSAMP 

10 

Y-HAT  -1.0217 
s\ 

Y-BAR 

-1.0000 

SIC-Y 

-0.0381 

Y-SHARP 

* 

0.0381 

y  given  x1 

plotted  In  Figure  6.7 
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APPENDIX  2.3: 


STATISTICS  FOR  HATCHED  FILTER  METHOD,  CENTER 
ELECTRODE,  NO  STIMULUS— FOR  MALE  SUBJECT 


5  R  REGMOD 
LU  OF  DATA  FILE 
90 

It  OF  INDEPENDENT  VARIABLES,  E.G.  3 
3 

It  OF  DIFFERENT  STIMULUS  CONDITIONS,  E.G.  2 
2 

0  OF  POINTS  IN  GROUP,  E.G.  10 
10 


TEMP1K.OUT 

1 

TEMP2K . OUT 

2 

TEMP3K. OUT 

3 

TEMP IK. OUT 

1 

TEMP 2K. OUT 

2 

TEMP3K.OUT 

3 

*1 

x2 

X3 

y 

0.1864E+01 

0.1881E+01 

0. 1748E+01 

0.1000E  *-01 

0. 1699E+01 

0. 2101E+01 

0. 2152E+01 

0. 100QE+01 

0. 1713E+01 

0. 1763E+01 

0. 1920E+01 

0. 1000E+01 

0.1965E+01 

0.2074E+01 

0. 1985E+01 

0 . 1000E+01 

0. 1672E+01 

0. 1302E+01 

0.1936E+01 

0.1000E+01 

0.2297E+01 

0.2983E+0I 

0. 1990E+01 

0.1000E+01 

0.1999E+01 

0.1824E+01 

0. 2178E+01 

0.1000E+01 

0. 1763E+01 

0.1483E+01 

0. 1953E+01 

0. 1000E+01 

0. 2073E+01 

0.2428E+01 

0. 1824E+01 

0. 1000E+01 

0.1558E+01 

0. 1715E+01 

0. 1882E+01 

0 . 1000E+01 

-0.1849E+01 

0.1645E+00 

0.2773E-01 

O.OOOOE+OO 

0. 2925E-01 

0. 2912E+00 

-0.25S9E+00 

O.OOOOE+OO 

0. 6134E-01 

-0.1876E-01 

-0. 2101E-01 

O.OOOOE+OO 

-0. 7798E-01 

0.2049E-01 

-0.7963E-01 

O.OOOOE+OO 

-0.8653E-01 

-0.6413E-01 

-0.7649E-01 

O.OOOOE+OO 

0.1701E+00 

-0. 9510E+00 

-02 . 633E+00 

O.OOOOE+OO 

0.1413E+00 

0.3667E+00 

0. 1326E+00 

O.OOOOE+OO 

0.I007E+00 

-0.119 5E+00 

-0.1432E+00 

O.OOOOE+OO 

0.2050E+00 

0.2188E+00 

0. 5784E-02 

O.OOOOE+OO 

-0. 6515E-01 

-0.1447E+00 

0. 2094E-01 

O.OOOOE+OO 

GIVE  STIMULUS  CONDITION,  E.G.  Y-l 


(verxjT) 


>  selected 


REG-CONSTANT  +  COEFFICIENTS 
-.02535  .55115  0  0 

neglected 

—  .  i 


Y-HAT 

* 

- 

■0.0254  + 

0.5512  XI  + 

0.0000  X2 

+  0. 

0000  X3  + 

S3  AMP 

=* 

1 

Y-HAT 

-0.0355 

Y-BAR 

* 

-0.0355SIC-Y 

-o.oooo 

Y -SHARP 

* 

0.0000 

NSAMP 

a 

2 

Y-HAT 

*- 

-0.0092 

Y-BAR 

-0.0224  SIG-Y 

-0.0186 

Y-SHARP 

* 

-0.8312 

N3A.'!P 

= 

3 

Y-HAT 

a 

0.0085 

Y-BAR 

*- 

-0. 0121 SIG-Y 

-0.0221 

Y-SHARP 

» 

-1.8292 

NSAMP 

= 

4 

Y-HAT 

— 

-0.0683 

Y-BAR 

-0.0262SIG-Y 

-0.0334 

Y-SHARP 

* 

-1.2777 

NSAMP 

= 

5 

Y-HAT 

—0.0730 

Y-BAR 

* 

-0.0355 SIG-Y 

•0.0357 

Y-SHARP 

- 

-1.0058 

NSAMP 

a 

6 

Y-HAT 

0.0684 

Y-BAR 

-0.0182 SIG-Y 

0.0531 

Y-SHARP 

-2.9174 

NSAMP 

a 

7 

Y-HAT 

0.0525 

Y-BAR 

-0.0081 SIG-Y 

■0.0554 

Y-SHARP 

- 

-6.8311 

NSAMP 

a 

B 

Y-HAT 

0.0302 

Y-BAR 

-0. 0033  SIG-Y 

-0.0530 

Y-SHARP 

- 

-15.9608 

NSAMP 

9 

Y-HAT 

0.0877 

Y-BAR 

- 

0. 0068SIC-Y 

-0.0581 

Y-SHARP 

- 

8 . 5662 

NSAMP 

= 

10 

Y-HAT 

-0.0613 

Y-BAR 

- 

0. 0000 SIG-Y 

-0.0589 

Y-SHARP 

* 

-3349.4158 

9  given  Xy 

plotted  in  figure  6,8 
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APPENDIX  2.4:  STATISTICS  FOR  MATCHED  FILTER  METHOD,  CENTER 
ELECTRODE,  WITH  STIMULUS — FOR  MALE  SUBJECT 


$  R  REGMOD 
LU  OF  DATA  FILE 
90 

t  OF  INDEPENDENT  VARIABLES  3 

3Y 

1  OF  DIFFERENT  STIMULUS  CONDITIONS  ,  E.G.  2 

2 

#  OF  POINTS  IN  GROUP,  E.C.  10 
10 

TEMP1K . OUT  1 

TEMP2K .  OUT  2 

TEMP3K.0UT  3 

TEMP1K.0UT  1 

TEMP 2K. OUT  2 

TEMP3K.0UT  3 


*1 

x2 

x3 

y 

0.1864E+01 

0. 1881E+01 

0. 1748E+01 

0 . 1 OOOE+Ol 

0. 1699E+01 

0. 2101E+01 

0.  2L52E+01 

0. 1000E+01 

0.1713E+01 

0. 1763E+01 

0. 1920E+01 

0. 1000E+01 

0.1965E+01 

0.2074E+01 

0. L98SE+01 

0.1000E+01 

0.1672E+01 

0.130XE+01 

0. L936E+01 

0. 1000E+01 

0.2297E+01 

0. 2983E+01 

0. 1990E+01 

0 . 1000E+01 

0. 1999E+01 

0.1824E+01 

O'.  2178E+01 

0. 1 OOOE+Ol 

0.1763E+01 

0.1483E+01 

0. 19S3E+01 

0.1 OOOE+Ol 

0.2073E+01 

0.2428E+01 

0. 18 24 Ml 

0. 1 OOOE+Ol 

0. 1558E+01 

0. 1715E+01 

0. 1882E+01 

0.1 OOOE+Ol 

-0.1849E+01 

0. 1645E+00 

0.27736-01 

0 . OOOOE+OO 

0. 2925E-01 

0.2912E+00 

-0. 2  559E+00 

0. OOOOE+OO 

0. 6134E-01 

-0. 1876E-01 

-0. '101E-01 

0. OOOOE+OO 

-0. 7798E-01 

0 , 2049E-01 

-0.  7 963E-01 

0. OOOOE+OO 

-0.8653E-01 

-0.6413E-01 

-0. 7649E-01 

0. OOOOE+OO 

0. 1701E+00 

-0.9510E+00 

-02 . 633E+00 

0. OOOOE+OO 

0.1413E+00 

0. 3667E+00 

0. 1326E+00 

0. OOOOE+OO 

0. 1007E+00 

-0.U95E+00 

-0. 1432E+00 

0.  OOOOE+OO 

0.2050E+00 

0.2188E+00 

0. 5784E-02 

0 . OOOOE+OO 

-0. 6515E-01 

-0.144 7 E+00 

0.  .'094E-01 

0. OOOOE+OO 

selected 


GIVE  STIMULUS  CONDITION,  E.G.  i-l 
1 

REG-CONSTANT  +  COEFFICIENTS 
-.02535  .55115  0  0 


neglected 


Y-HAT 

» 

-0 

.0254  + 

0.5512  :.i  + 

0.0000  X2 

+  0. 

0000  X3  + 

N3AMP 

- 

1 

Y-HAT 

1.0023 

Y  BAR 

1.0023 

SIG-Y 

-  0.0000 

Y-SHARP 

-  0.0000 

NSAMP 

2 

Y-HAT 

0.9110 

Y-BAR 

0.9566 

SIG-Y 

*  0.0646 

Y-SHARP 

*0.0675 

N3AMP 

* 

3 

Y-HAT 

0.9188 

Y-BAR 

0.9440 

SIG-Y 

-  0.0506 

Y-SHARP 

-0.0536 

NSAMP 

« 

4 

Y-HAT 

1.0576 

Y-BAR 

0.9724 

SIG-Y 

-0.0702 

Y-SHARP 

-  0.0722 

NSAMP 

* 

5 

Y-HAT 

0.8964 

Y-BAR 

0.9572 

SIG-Y 

-0.0697 

Y-SHARP 

-0.0728 

NSAMP 

6 

Y-HAT 

1.2405 

Y-BAR. 

1.0044 

SIG-Y 

-0.1314 

Y-SHARP 

-0.1308 

NSAMP 

-= 

7 

Y-HAT 

1.0766 

Y-BAR 

1.0147 

SIG-Y 

-  0.1230 

Y-SHARP 

-0.1212 

NSAMP 

- 

8 

Y-HAT 

0.9463 

Y-BAR 

1.0062 

SIG-Y 

-0.1164 

Y-SHARP 

-0.1157 

NSAMP 

* 

9 

Y-UAT 

1.1173 

Y-BAR 

1.0185 

SIG-Y 

-  0.1150 

Y-SHARP 

-0.1129 

NSAMP 

* 

10 

Y-HAT 

0.8335 

Y-BAR 

1.0000 

SIC-Y 

-  0.1232 

Y-SHARP 

-0.1232 

A 

y  given  Xj 

plotted  In  Figure  6.8 


model 

values 


APPENDIX:  3.1:  STATISTICS  FOR  MATCHED  FILTERS,  ALL  ELECTRODES, 
BOTH  STIMULI —FOR  FEMALE  SUBJECT 


TEMPI 0.001 

1 

TEMP2D.0UT 

2 

TEMP 3D. OUT 

3 

TEMP1D.0UT 

1 

TEMP2D.0UT 

2 

TEMP3D.0UT 

3 

*1 

*2 

X3 

0. 9931E-03 

0. 6960E-01 

0.3915E-01 

0.3194E-01 

0. 6296E-01 

0. 1829E-02 

0.6406E-01 

-0.8106E-02 

-0.187  5E+00 

0. 2677E-02 

0.2856E-01 

0.2  L38E+00 

0. 3547E-01 

0. 1780E-01 

-0. 1254E+00 

0 . 3902E-01 

-0. 3840E-01 

-0.4314E-01 

0. 1096E+00 

0.6040E-01 

-0. 84 J2E-01 

0.4794E-01 

0. 3290E-01 

-0. 6859E-01 

0.4677E-01 

0. 1066E+00 

0. 1987E+00 

0. 1385E-01 

0. 5961E-01 

0.  9772E-01 

0. 1 941E+01 

0.2001E+01 

0. 1639E+01 

0. 2072E+01 

0. 2103E+01 

0. 1995E+01 

0. 1939E+01 

0.2029E+01 

0. 1860E+01 

0. 2014E+01 

0.2056E+01 

0. 1869E+01 

0. 2068E+0L 

0.2039E+01 

0. 1 986E+01 

0.1826E+01 

0. 1993E+01 

0. L806E+01 

0. 2012E+01 

0. 1916E+01 

0.167 1E+01 

0. 1987E+01 

0. 1933E+01 

0.  2117E+01 

0.1934E+01 

0.1965E+01 

0. 2377E+01 

0 . 2024E+01 

0. 2081E+01 

0.2496E+01 

GIVE  STIMULUS  CONDITION 

.  E.C.  Y«1 

*-  o 


y 

O. OOOOE+OO 
0.  OOOOE+OO 
0. OOOOE+OO 
0. OOOOE+OO 
0. OOOOE+OO 
0. OOOOE+OO 
0. OOOOE+OO 
0. OOOOE+OO 
0.  OOOOE+OO 
0.  OOOOE+OO 
0.1000E+01 
0.  I OOOE+Ol 
0. lOOOE+Ol 
0. I OOOE+Ol 
0. 1 OOOE+Ol 
0. LOOOE+Ol 
0. 1 OOOE+Ol 
0. 1 OOOE+Ol 
0, lOOOE+Ol 
0.1 OOOE+Ol 


VERDJ123 


REG-CONSTANT  +  COEFFICIENTS 
.-.0198  .1736  .3332  -.0078 

-0.0198  •* 


Y-HAT 


0.1736  XI  + 


constant  response  In 
Yle£t  hemisphere- important 

*  n 


0.3332  X2  + 


-0.0078  X3  + 


MS AMP  -  L 

Y-HAT  -  0.0033 

Y-BAR 

*  0.0033  SIG-Y 

-0.0000 

Y-SHARP 

0.0000 

NSAMP  ■  2 

Y-HAT  -  0.0067 

Y-BAR 

-  0.0050  SIG-Y 

-0.0024 

Y-SHARP 

0.4894 

N3AMP  -  3 

Y-HAT  -  -0.0099 

Y-BAR 

*  0.0000  SIG-Y 

-0.0088 

Y-SHARP 

544.2001 

NSAMP  -  4 

Y-HAT  «  -0.0117 

Y-BAR 

* -0.0029  SIG-Y 

-  0.0093 

Y-SHARP 

-3.1736 

NSAMP  -  5 

Y-HAT  -  -0.0067 

Y-BAR 

-  -0.0037  SIG-Y 

-  0.0082 

Y-SHARP 

-2.2275 

NSAMP  -  6 

Y-HAT  -  -0.0255 

Y-BAR 

- -0.0073  SIC-Y 

*0.0115 

Y-SHARP 

-1.5767 

NSAMP  -  7 

Y-HAT  »  0.0200 

Y-BAR 

-  -0.0034  SIG-Y 

-  0.0147 

Y-SHARP 

-4.3235 

NSAMP  -  8 

Y-HAT  -  0.0000 

Y-BAR 

-  -0.0030  SIG-Y 

-0.0137 

Y-SHARP 

-4.5969 

NSAMP  -  9 

Y-HAT  -  0.0223 

Y-BAR 

* -0.0002  SIG-Y 

-0.0153 

Y-SHARP 

-87.7774 

IIS  AMP  -  10 

Y-HAT  -  0.0017 

Y-BAR 

-  0.0000 

SIG-Y 

-0.0145 

Y-SHARP 

1107.8656 

TERMINATE  * 

n 

-1,  SAME  FILE  - 

0,  NEW 

FILE  -  1 

GIVE  STIMULUS  CONDITION,  E 

C.  Y-l 

y  given  x 

y 

NSAMP  -  1 

Y-HAT  -0.9711 

Y-BAR 

-0.9711 

SIG-Y 

-0.0000 

Y-SHARP 

NSAMP  -  2 

Y-HAT  -1.0252 

Y-BAR 

-0.9981 

SIG-Y 

-0.0383 

Y-SHARP 

m 

0.0384 

NSAMP  *  3 

Y-HAT  -0.9782 

Y-BAR 

-0.9915 

SIG-Y 

-0.0294 

Y-SHARP 

m 

0.0297 

NSAMP  -  4 

Y-HAT  -1.0001 

Y-BAR 

-0.9937 

SIG-Y 

-0.0244 

Y-SHARP 

m 

0.0246 

NSAMP  -  5 

Y-HAT  -1.0033 

Y-BAR 

-0.9956 

SIG-Y 

-0.0216 

Y-SHARP 

*k 

0.0217 

NSAMP  -  6 

Y-HAT  -0.9472 

Y-BAR 

-0.9875 

SIC-Y 

-0.0276 

Y-SHARP 

■ 

0.0280 

(0 

a 

NSAMP  -  7 

Y-HAT  -0.9547 

Y-BAR 

•0.9828 

SIC-Y 

-0.0281 

Y-SHARP 

• 

0.0286 

a 

NSAMP  -  8 

Y-HAT  -0.9528 

Y-BAR 

-0.9791 

SIG-Y 

-0.0281 

Y-SHARP 

m 

0.0287 

I 

NSAMP  -  9 

Y-HAT  -0.9522 

Y-BAR 

-0.9761 

SIC-Y 

-0.0278 

Y-SHARP 

» 

0.0284 

4J 

(A 

NSAMP  -  10 

Y-HAT  -1.0055 

Y-BAR 

-0.9790 

SIG-Y 

-0.0278 

Y-SHARP 

* 

0.028 4 

TERMINATE  * 
-1 

-l,  SAME  FILE  - 

0,  NEW 

FILE  -  1 

X-SC 

Y-SC 

O.K. 

0 


000000 

000000 


5 . 000000 
1.000000 
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no  stimulus 


APPENDIX  3.2:  EXAMPLE  OF  AUGMENTING  DATA  FOR  APPROXIMATE 
LIKELIHOOD  ESTIMATION,  ALL  CHANNELS— FOR 
FEMALE  SUBJECT 


ASS  VERDJR.DAT  FOR090 

PREVIOUS  LOGICAL  NAME  ASSIGNMENT  REPLACED 
$  R  REGRESS 
LU  OP  DATA  FILE 
90 

#  OF  INDEPENDENT  VARIABLES  ,E.C.  3 
3 

II  OF  DIFFERENT  STIMULUS  CONDITIONS,  E.C  2 
3 


II  OF  POINTS  IN  GROUP,  E.G. 
10 

TEMP1D.OUT 

TEMP2D.OUT 

TEMP3D.OUT 

EMPOO 

EMPOO 

EMPOO 

EMPOO 

EMPOO 

EMPOO 


0.1941E+01  0.2001E+01 

0. 2072E+01  0.2103E+01 

0.1939E+01  0.2029E+01 

0.2014E+01  0.2056E+01 

0.2068E+01  0. 2039E+01 

0.1826E+OL  0.I993E+O1 

0. 2012E+01  0.19I6E+O1 

0.1987E+01  0. 1933E+01 

0. 1934E+01  0. 1965E+01 

0.2024E+01  0.2081E+01 

0 . OOOOE+OO  0 . OOOOE+OO 
0. OOOOE+OO  0. OOOOE+OO 
0. OOOOE+OO  0. OOOOE+OO 
0. OOOOE+OO  0. OOOOE+OO 
0 . OOOOE+OO  0 . OOOOE+OO 
0. OOOOE+OO  0. OOOOE+OO 
0 . OOOOE+OO  0 . OOOOE+OO 
0 . OOOOE+OO  0 . OOOOE+OO 
0 . OOOOE+OO  0. OOOOE+OO 
0. OOOOE+OO  0. OOOOE+OO 
0. OOOOE+OO  0. OOOOE+OO 
0 . OOOOE+OO  0 . OOOOE+OO 
0 . OOOOE+OO  0 . OOOOE+OO 
0. OOOOE+OO  0. OOOOE+OO 
0. OOOOE+OO  0. OOOOE+OO 
0. OOOOE+OO  0. OOOOE+OO 
0. OOOOE+OO  0. OOOOE+OO 
0. OOOOE+OO  0. OOOOE+OO 
0. OOOOE+OO  0. OOOOE+OO 
0. OOOOE+OO  0. OOOOE+OO 
SIGNIFICANCE  ALFA,  E.G..05 
.05 

BECOVM  IER  *  0 

A-MATRIX 
0. 2623E+02 

0. 2659E+02  0.2701E+02 

0. 2622E+02  0.2661E+02 

0.1321E+02  0.1341E+02 


10 

1 

2 

3 

0 

0 

0 

0 

0 

0 


0. 1639E+01  0. 1000E+Q1 
0.1995E+01  0.1000E+01 
0.18  60E+01  0 . 1000E+01 
0. 1869E+01  0. 1000E+01 
0. 1986E+01  Q.1000E+01 
0. 1806E+01  0.1000E+01 
0.1671E+01  0. 1000E+01 
0.2117E+01  0. 1000E+01 
0.2377E+01  0.1000E+01 
0.2496E+01  0. 1000E+01 
0 . OOOOE+OO  0 . OOOOE+OO 
0. OOOOE+OO  0. OOOOE+OO 
0. OOOOE+OO  0. OOOOE+OO 
0. OOOOE+OO  0. OOOOE+OO 
0. OOOOE+OO  0. OOOOE+OO 
0. OOOOE+OO  0. OOOOE+OO 
0. OOOOE+OO  0. OOOOE+OO 
0. OOOOE+OO  0. OOOOE+OO 
0 . OOOOE+OO  0 . OOOOE+OO 
0. OOOOE+OO  0. OOOOE+OO 
0. OOCOE+OO  0. OOOOE+OO 
0. OOOOE+OO  0. OOOOE+OO 
0. OOOOE+OO  0. OOOOE+OO 
0. OOOOE+OO  0. OOOOE+OO 
0. OOOOE+OO  0. OOOOE+OO 
0. OOOOE+OO  0. OOOOE+OO 
0. OOOOE+OO  0. OOOOE+OO 
0. OOOOE+OO  0. OOOOE+OO 
0. OOOOE+OO  0. OOOOE+OO 
0. OOOOE+OO  0. OOOOE+OO 


0. 2689E+02 

0. 1321E+02  0. 6667E+01 


regress  DJ 


added  points  at  origin  In  order 
to  approximate  maximum  likelihood 
approach 
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(Cont’d  on  next  page) 


APPENDIX  3.2  (Cont’d.) 


RLMUL  ier  »  0 

DF  DF-RES  DF-TOT  SS-REG  SS-RES  SS-TOT  MS-REG 

3.0000  26.0000  29.0000  0.6660E+01  0.6952E-02  0.6667E+01  0.2220E+01 

MS-RES  F-VALUE  P-EXC  HO  X-EXPL  STDV-RES  Z-R-MEAN  0K-4 


0. 2674E-03  0.8303E+04 


0.163SE-01  0.4906E+01 


B-MATRIX 

REG-COEFF 

0. 1736E+00 
0. 3332E+00 
-0. 7805E-02 
0.3476E-03 


STDV-ERR 


PART- F-VALUE 


0. 2250E-01  0.3248E+00  0.7353E-01  0.1491E-02 
0.1829E+00  0.4836E+00  0.7315E-01  0.5548E-02 
-0.4906E-01  0.3345E-01  0.2007E-01  0.4044E-04 
-0. 7167E-02  0.7862E-02  0.3656E-02  O.OOOOE+OO 


0. 5577E+01 
0.207 5E+02 
0.1512E+00 
O.OOOOE+OO 


VARB 

0. 2022E+02 

-0. 1936E+02  0.2002E+02 

-0. 5575E+00  -0. 9352E+00  0.1506E+01 


0.1736  XI  + 


may  be  slightly 
readjusted 


0.3332  X2  + 


model  values  used 


model  DJ 


-0.0078  X3  + 
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APPENDIX  3.3:  STATISTICS  FOR  MATCHED  FILTERS, 
BOTH  STIMULI— FOR  MALE  SUBJECT 


TEMP1K.OUT 

1 

TEMP2K.OUT 

2 

TEMP3K.OUT 

3 

TEMP1K.OUT 

1 

TEMP 2K. OUT 

2 

TEMP3K.0UT 

3 

*1 

X2 

x3 

y 

0.1864E-t-01 

0. 1881E+01 

0. 1748E+01 

0.1000E+01 

0. 1699E+01 

0.2101E+01 

0.2152E+01 

0.1000E+01 

0. 1713E+01 

0. 1763E+01 

0. 1920E+01 

0. 1000E+01 

0. 1965E+01 

0.2074E+01 

0. 1985E+01 

0 . 1000E+01 

0. 1672E+01 

0. 1302E+01 

0. 1936E+01 

0.1000E+01 

0. 2297E+01 

0. 2983E+01 

0.1990E+01 

0.1000E+01 

0. 1999E+01 

0. 1824E+01 

0.2178E+01 

0. 1000E+01 

0. 1763E+01 

0. 1483E+01 

0.1953E+01 

0. 1000E+01 

0.2073E+01 

0. 2428E+01 

0. 1824E+01 

0. 1000E+01 

0. 1558E+01 

0. 1715E+01 

0. 1882E+01 

0.1000E+01 

-0. 1849E+01 

0. 1645E+00 

0.2773E-01 

0. OOOOE+OO 

0. 2925E-01 

0. 2912E+00 

-0. 2559E+00 

0. OOOOE+OO 

0. 6134E-01 

-0. 1876E-01 

-0. 2101E-01 

0. 0000E+00 

-0. 7798E-01 

0. 2049E-01 

-0.7  963E-01 

0. OOOOE+OO 

-0. 8653E-01 

-0.6413E-01 

-0. 7649E-01 

0. OOOOE+OO 

0. 1701E+00 

-0.9510E+00 

-02 . 633E+00 

0. OOOOE+OO 

0. —13E+00 

0. 3667E+00 

0. 1326K+00 

0. OOOOE+OO 

0.1007E+00 

-0.  U95E+00 

-0.1432EH30 

0. OOOOE+OO 

0.2050E+00 

0.2188E+00 

0.  5784E-02 

0 . OOOOE+OO 

-0. 6515E-01 

-0.144 7 E+00 

0.2094E-01 

0. OOOOE+OO 

GIVE  STIMULUS 

CONDITION, 

E.G.  Y-l 

RRG-CONSTANT  +  COEFFICIENTS 
.0X33  .2077  -.0550  .3662 


/^Y-HAT  - 

0.0133 

+  0.2077 

XI  +  - 

■0.0550  X2  +  0. 

3662  X3 

N3AMP 

3 

1 

Y-HAT 

-  0.0106 

Y-BAR 

-  0.0l06SlG-Y 

-  0.0000 

model 

NSAMP 

2 

Y-HAT 

-  -0.0903 

Y-BAR 

-  -0.0399 SIG-Y 

-  0.0714 

values 

NSAMP 

3 

Y-HAT 

-  0.0194 

Y-BAR 

-  -0.0201 SIG-Y 

»  0.0610 

NSAMP 

4 

Y-HAT 

-  -0.0359 

Y-BAR 

-  -0.0241  SIG-Y 

-  0.0504 

NSAMP 

5 

Y-HAT 

-  -0.0292 

Y-BAR 

-  -0.0251  SIG-Y 

-  0.0437 

NSAMP 

6 

Y-HAT 

-  0.0045 

Y-BAR 

-  -0.0202  SIG-Y 

-  0.0409 

NSAMP 

7 

Y-HAT 

-  0.0710 

Y-BAR 

-  -0.0071 SIG-Y 

-  0.0508 

NSAMP 

8 

Y-HAT 

-  -0.0116 

Y-BAR 

-  -0.0077  SIG-Y 

>  0.0471 

NSAMP 

9 

Y-HAT 

-  0.0460 

Y-BAR 

-  -0.0017  SIG-Y 

-  0.0475 

NSAMP  -  10  Y-HAT  -  0.0154 
TERMINATE  -  -1,  SAME  FILE  -  ( 
0 

GIVE  STIMULUS  CONDITION,  E.G 
1 

Y-BAR 
3,  NEW 

.  Y-l 

-  0.0000  SIG-Y 
FILE  -  1 

-  0.0451 

NSAMP 

l 

T-HAT 

-0.9373 

Y-BAR 

=■0.9373  SIG-Y 

-0.0000 

NSAMP 

2 

Y-HAT 

-1.0385 

Y-BAR 

*0.9879  SIG-Y 

-0.0716 

NSAMP 

3 

Y-HAT 

-0.9753 

Y-BAR 

*0.9837  SIG-Y 

-0.0512 

NSAMP 

4 

Y-HAT 

-1.0341 

Y-BAR 

-0.9963  SIG-Y 

-0.0488 

NSAMP 

5 

Y-HAT 

-0.9978 

Y-BAR 

*0.9966  SIG-Y 

-0.0423 

NSAMP 

6 

Y-HAT 

-1.0549 

Y-BAR 

-1 . 0063  SIG-Y 

-0.0447 

NSAMP 

7 

Y-HAT 

-1.1260 

Y-BAR 

-1.0234  SIG-Y 

•0.0609 

NSAMP 

8 

Y-HAT 

-1.0132 

Y-BAR 

-1.0221  SIG-Y 

-0.0565 

NSAMP 

9 

Y-HAT 

-0.9783 

Y-BAR 

-1.0173  SIG-Y 

-0.0548 

NSAMP  -  10 
TERMINATE  - 
-1 

Y-HAT  -0.9318 
-1,  SAME  FILE  - 

Y-BAR 
0,  NEW 

-1.0087  SIC-Y 
PILE  -  1 

-0.0583 

X-SC  1.000000  5.000000 

Y-SC  -1.000000  1.000000 

O.K.-l 
0 


ALL  ELECTRODES, 


VERKJ123 


constant  response  In  right 
hemisphere  -  important 

+  N, 

Y-SHARP  -  0.0000 
Y-SHARP  - -1.7889 
Y-SHARP  - -3.0281 
Y-SHARP  - -2.0928 
Y-SHARP  - -1.7415 
Y-SHARP  - -2.0300 
Y-SHARP  - -7.1292 
Y-SHARP  - -6.1200 
Y-SHARP  - -27.4603 
Y-SHARP  - -2442.5129 


Y-SHARP  - 
Y-SHARP  -  0.0725 
Y-SHARP  -  0.0520 
Y-SHARP  -  0.0490 
Y-SHARP  -  0.0424 
Y-SHARP  «  0.0444 
Y-SHARP  -  0.0595 
Y-SHARP  -  0.0553 
Y-SHARP  -  0.0539 
Y-SHARP  -  0.0578  .. 


stimulus  no  stimulus 


APPENDIX  3.4:  EXAMPLE  OF  AUGMENTING  DATA  FOR  APPROXIMATE  MAXIMUM 

LIKELIHOOD  ESTIMATION.  ALL  CHANNELS — FOR  MALE  SUBJECT 


OF  INDEPENDENT  VARIABLES,  E.G.  3 
OF  DIFFERENT  STIMULUS  CONDITIONS,  E.G.  2 


a  OF  POINTS  IN  GROUP*  E.G 
10 

.  10 

regress 

TEMP1K.0UT 

1 

TEMP2K.0UT 

2 

TEMP 3 K. OUT 

3 

EMPOO 

0 

EMPOO 

0 

EMPOO 

0 

EMPOO 

0 

EMPOO 

0 

EMPOO 

0 

0, 1864E+01 

0.188IE+01 

0. 1748E+01 

0. 1000E+01 

0.1699E+01 

0.2101E+01 

0. 2152E+01 

0. 1000E+01 

0.1713E+01 

0. 1763E+01 

0. 1920E+01 

0. 1000E+01 

0. 1965E+01 

0. 2074E+01 

0. 1985E+01 

0. 1000E+01 

0. 1672E+01 

0. I302E+01 

0. 1936E+01 

0. 1000E+01 

0. 2297E+01 

0. 2983E+01 

0. 1990E+01 

0. 1000E+01 

0.1999E+OI 

0. 1824E+01 

0. 2178E+01 

0. 1000E+01 

0. 1763E+0I 

0. 1483E+01 

0. I953E+0I 

0. lOOOE+Ol 

0.2073E+01 

0.  2428E+01 

0. 1824E+01 

0. 1000E+01 

0.1558E+01 

0. 17I5E+01 

0. 1882E+01 

0. 1000E+01 

O.OOOOE+OO 

0 . OOOOE+OO 

0. OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

0. OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

0 . OOOOE+OO 

O.OOOOE+OO 

0 . OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

SIGNIFICANCE 

ALFA,  E.G. . 

05 

.05 

BECOVM  IER  = 

■  0 

A-MATRIX 

0. 2352E+02 

0. 2504E+02 

0. 2757E+02 

0.2429E+02 

0. 2553E+02 

0. 2568E+02 

0.1240E+02 

0. 1303E+02 

0. 1304E+02 

0.666 7 E+01 

RLMUL  IER  - 

0 

ANOVA 

DF 

DF-RES 

DF-TOT  SS- 

REG  SS' 

added  points  at  origin 


3 . 0000 

MS-RES 
0. 1I83E-02 


26.0000 


29.0000  0. 6636E+01  0.3077E-0I 


SS-TOT 

0.6667E+01 


MS-REG 
0. 2212E+01 


F-VALUE 

0.1869E+04 


P-EXC  HO  Z-EXPL  STDV-RES  Z-R-MEAN  OK-4 

0.0000  99.5385  0.3440E-0I  0.1032E+02  4.0000 


(Cont'd  on  next  page) 
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APPENDIX  3.4  (Cont’d.) 


B-MATRIX 


REG-COEFF 
0.2077E+00 
-0. 5504E-01 
0.3662E+00 
0. 1539E-02 
VARB 

0. 5251E+01 
-0.2145E+01 
-0.2834E+01 


L-LIM 

0.4567E-01 

-0.1367E+00 

0.2657E+00 

-0.1426E-01 


0.1335E+01 
0. 7014E+00 


U-LIM 
0.3698E+00 
0. 2665E-01 
0.4668E+00 
0.1734E-01 


0.2022E+01 


STDV-ERR 
0. 7683E-01 
0. 3974E-01 
0.4891E-01 
0.7686E-02 


PART- F- VALUE 
0.8216E-02 
0.2270E-02 
0. 6634E-01 
0.  OOOOE+OO 


P-EXC  HO 
0.6942E+01 
0. 1918E+01 
0. 5605E+02 
0. OOOOE+OO 


Y-HAT  -  0.0015  +  0.2077  XI  +  -0.0550  X2  +  0.3662  X3  + 


J 


model 


model  values  used 


I 


APPENDIX  4.1s  USEFULNESS  OF  VARIABILITY  MEASURE — 
FOR  FEMALE  SUBJECT 


R  REGRESS 
LU  OF  DATA  FILE 
90 

If  OF  INDEPENDENT  VARIABLES,  E.G.  3 


If  OF  DIFFERENT  STIMULUS  CONDITIONS, 
2 


E.G,  2 


»  OF  POINTS  IN 

GROUP,  E.G. 

IQ 

10 

TEMPI D. OUT 

1 

TEMP2K.OUT 

2 

TEMP3D.OUT 

3 

TEMP  ID.  OUT 

1 

TEMP2K.OUT 

2 

TEMP3D.OUT 

3 

B, 

Z1 

Z2 

z3 

V 

0.3747E-02 

0.3559E-02 

0.1104E-01 

0. OOOOE+OO 

0. 3873E-02 

0. 1041E-01 

0.4610E-02 

0. OOOOE+OO 

0. 1695E-02 

0. 5232E-02 

0. 1454E-01 

0. OOOOE+OO 

0.8106E-02 

0. 5737E-02 

0. 1219E-01 

0. OOOOE+OO 

0. 7576E-02 

0.4068E-02 

0.1013E+00 

0. OOOOE+OO 

0. 9514E-02 

0. 5407E-02 

0. 1111E-01 

0. OOOOE+OO 

0. 1090E-01 

0. 6532E-02 

0.1675E-01 

0. OOOOE+OO 

0.1369E-01 

0. 1804E-02 

0.9163E-02 

0. OOOOE+OO 

0.1902E-02 

0.9029E-02 

0. 9211E-01 

0. OOOOE+OO 

0. 7406E-02 

0.  5252E-02 

0.2258E+00 

0. OOOOE+OO 

0. 6244E-02 

0.1617E-01 

0.2427E-01 

0.1000E+01 

0. 4223E-01 

0. 2207E-01 

0.2985E-01 

0. 1000E+01 

0. 1353E-01 

0.6600E-01 

0.2260E+00 

0. 1000E+01 

0.3528E-01 

0. 1180E-01 

O.1693E-01 

0. 1000E+01 

0. 5051E-01 

0.6275E-02 

0.2333E+00 

0.1000E+01 

0.3528E-01 

0. 1105E-01 

0.4538E-01 

0. 1000E+01 

0.1457E-01 

0.2269E-01 

0.5961E-01 

0. 1000E+01 

0. 1958E-01 

0. 1376E-01 

0.1881E+00 

0. 1000E+01 

0. 2853E-01 

0. 4474E-01 

0. 6369E-01 

0. 1000E+01 

0.2771E-01 

0.1310E-03 

0.5115E-01 

0, 1 OOOE+Ol 

VERDJW 


stimulus  condition 


.05 

BECOVM  IER 
A-MATRIX 
0.3990E-02 
0.6276E-03 
0.4885E-02 
0. 1025E+00 
RLMUL  IER  - 
ANOVA 

DF  DF-RES 

3.0000  16.0000 


0.4802E-02 
0.82 57E-02 
0. 7883E-01 


DF-TOT 

19.0000 


0.1224E+00 
0. 2198E+00 
0 

SS-REG 

0.3471E+01 


z  are  the  10-sec  segment 
variability  measures 


0.5000E+01 


SS-RES 
0. 1529E+01 


SS-TOT 
0 . 5000E+01 


MS-REG 

0.1157E+01 


MS- RES 
0.9557E-01 

B-MATRIX 


F-VALUE  P-EXC  HO/  X-EXPL)  STDV-RES  X-R-KEAN  OK-4 
0. 1211E+02  0.0002  (69.4173  0.3901E+00  0.6183E+02  4.0000 

— z_variability  z_  is  related  to  stimulus 


0. 2365E+02 
0. I342E+02 
-0. 5328E-01 
-0.8277E-01 
VARB 

0.2650E+03 
-0. 1861E+02 
-0-9320E+01 


0.1298E+02 
0.3330E+01 
-0. 2080E+01 
-0.3443E+00 


0.2369E+03 

-0.1524E+02 


0.3432E+02 
0. 2350E+02 
0, 1974E+01 
0, 1788E+00 


0. 9568E+01 


0.5032E+01 
0.4758E+01 
0. 9563E+00 
0.1 234E+00 


^13.4174  X2 


012111E+01 
0.7599E+00 
0. 2966E-03 
0 . OOOOE+OO 


D.2209E+02 
0. 7952E+01 
0.3104E-02 
0. OOOOE+OO 


Y-HAT  -  -0.0828  + 


23.6526  XI  + 
center 


left  hemisphere  right 
more  variable 
than  right 


+  -0.0533  X3  + 

) 
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APPENDIX  4.2:  USEFULNESS  OF  VARIABILITY  MEASURE— 
FOR  MALE  SUBJECT 


LU  OF  DATA  FILE 
90 

H  OF  INDEPENDENT  VARIABLES,  E.G.  3 
3 

II  OF  DIFFERENT  STIMULUS  CONDITIONS  .  E.G.  2 
2 

II  OF  POINTS  IN  GROUP  ,  E.G.  10 
10 

TEMP1K.0UT 


TEMP2K.0UT 
TEMP3K.0UT 
TEMP1K.0UT 
TEMP2K.0UT 
TEMP3K. OUT 


1 

Z2 

3 

y 

0.9354E+00 

0.5004E-01 

0. 1019E+00 

0. OOOOE+OO 

0.7964E+00 

0.1912E+00 

0. 6059E-01 

0.0000  00 

0.7046E+00 

0.3850E+00 

0.4898E-01 

0. OOOOE+OO 

0.1252E+00 

0.2421E+00 

0. 1624E-0. 

0. OOOOE+OO 

0. 1737E+01 

0.1029E+00 

0.4250E-01 

0. OOOOE+OO 

0.3046E+00 

0.1029E+01 

0. 5255E-01 

0. OOOOE+OO 

0.4322E+00 

0. 1029E+01 

0.4749E- 01 

0. OOOOE+OO 

0.8312E+00 

0.3  3E+00 

0.2058E+00 

0. OOOOE+OO 

0.8436E+00 

0.2217E+00 

0.2985E-01 

0. OOOOE+OO 

0.1105E+01 

0.2338E-01 

0. 1348E+O0 

0. OOOOE+OO 

0. 5642E-01 

0. 3131E+00 

0. 3875E-01 

O.IOOOE+Ol 

0.1492E+00 

0. 1244E+00 

0. 1989E-01 

O.IOOOE+Ol 

0. 1652E+00 

0. 6974E+00 

0.4122E-01 

O.IOOOE+Ol 

0.1554E+00 

0.8784E-01 

0.1984E-01 

O.IOOOE+Ol 

0. 9855E+00 

0.1420E+00 

0.2439E+OO 

0. 1000E+01 

0. 5516E+00 

0. 1642E+01 

0.1220E+00 

0. 1000E+01 

0. 3278E-01 

0.1367E+01 

0.2906E+00 

0. 1000E+01 

0. 3021E+00 

0.2310E+00 

0. 1102E+00 

0. 1000E+01 

0.  7830E-00 

0.8025E+00 

0.7498E-01 

O.IOOOE+Ol 

0. 5214E+00 

0.8611E+00 

0. 6018E-01 

O.IOOOE+Ol 

.05 

BECOVM  IER  - 
A-MATRIX 
0.3889E+01 
-0. 1068E+00 
0. 1428E+00 
-0. 2444E+01 
RLMUL  IER  - 
ANOVA 

DF  DF-RES 
3.0000  16.0000 


0.3975E+01 
0.2379E+00 
0. 1488E+01 


0.8428E-01 
0. 7495E-01 


0. 5000E+01 


VERKJW 


DF-TOT 

19.0000 


SS-REG 
0. 2139E+01 


MS-RES 
0. 1788E+00 


F-VALUE 
0. 3988E+01 


P-EXC  HO/  %-EXPL 
0.0268  (42.7825. 


SS-RES 

0.2861E+01 

STDV-RES 
0.4229E+00 


SS-TOT 
0 . 5000E+01 

X-R-MEAN 
0.8457E+02 


MS-REG 

0.7130E+00 


OK-4 

4.0000 


B-MATRIX 

REG-COEFF  L-LIM  U-LIM 

-0. 6653E+00  -0. 1140E+01  -0.1908E+00 
0. 2836E+00  -0.2150E+00  0.7822E+00 
0.1216E401  -0.2219E+01  0.4750E+01 

0. 6175E+00  0. 2108E+00  0.1024E+01 
VARB 

0 . 2803E+00 

0.4226E+01  0. 3094E+00 

-0. 5969E+00  -0. 9468E+00  0.1555E+02 

Y-HAT  -  0.6175  +  -0.6653  XI  + 

center 


variability  z  Is  related  to  stlaulus 


STDV-ERR 
0.2339E+00 
0.2352E+00 
0.1667E+01 
0. 1918E+00 


0.2836  X2  + 
left 


PART-F-VALUE  P-EXC  HO 


0.1580E+01 
0.2599E+00 
0.9506E-01 
0 . OOOOE+OO 


0.8834E401 
0..1454E-01 
0.5316E400 
0. OOOOE+OO 


J..2157  X3,+ 

right  healsphere  taore  variable  than  left 
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